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UNDER  DYNAMIC  AND  QUASI-STATIC  CYCLIC  LOADING 

by 

A.M.  Eleiche 


ABSTRACT 


In  Part  I  of  the  present  investigation,  the  modified  version  of  the 
torsional  split  Hopkinr/on  bar  apparatus  was  further  exploited  to  generate 
new  data  bases  for  OFHC  copper  and  mild  steel  related  to  their  stress- 
strain  characteristics  at  ambient  temperature  under  sequential  reverse 
loading  at  impact  strain  rates,  and  different  large  values  of  the  constant 
cyclic  strain  amplitude.  The  obn^rved  response  was  compared  with  that 
obtained  under  siviilar  but  quasi-static* strain  rate  conditions.  All  data 
were  evaluated  in  terms  of  various  Bauscbinger  effect  parameters.  It  is 
shown  that  the  common  practice  of  using  quasi-static  reverse  loading  charac¬ 
teristics  in  the  analysis  of  problems  where  impact  strain  rates  are  predomin¬ 
ant  can  yield  to  an  overly  safe  and  non-economical  design,  particularly  if 
strain  rate  sensitive  materials  are  being  used.  As  expected,  the  present 
data  do  not  conform  with  the  predictions  of  classical  plasticity  theories, 
since  experimental  conditions  involved  essentially  large,  high  strain-rate 
and  adiabatic  deformations.  An  examination  of  the  Currently  widening litera¬ 
ture  on  advanced  constitutive  equations  for  metals ^through  many  of  the  recent 
comparative  reviews  seems  to  indicate  that  the  theories  of  Bodner  and 
associates  and  of  Cernocky  and  Krerapl  offer  the  greatest  promises  for  model¬ 
ing  the  present  data.  Further  in-depth  studies  in  that  direction  are  being 
recommended.  ' 

In  Parc  II,  the  problem  of  the  simple  thick-walled  tube  subjected  to 
steady  internal  pressure  and  cyclic  thermal  gradient,  was  investigated.  A 
new  numerical  procedure  was  constructed  to  solve  thesame  problem.  The 
validity  of  this  procedure  was  extensively  verified  through  comparisons 
with  available  analytical  solutions  based  on  the  kinematic  hardening  assump¬ 
tion  for  an  elastic-linearly  hardening  tube  material  and  the  uniaxial  model 


of  the  tube.  This  numerical  technique  is  further  elaborated  to  accomodate 
the  isotropic  hardening  rule  and  to  utilize  a  more  generalized  biaxial  tube 
model.  This  model  is  found  to  be  more  accurate  and  superior  to  the  uniaxial 
model  being  widely  used  in  similar  analyses. 


FOREWORD 


This  technical  report  was  prepared  by  Professor  A.M.  Eleiche  of  the 
Department  of  Mechanical  Design  and  Production,  Faculty  of  Engineering, 

Cairo  University, Egypt.  The  work  reported  herein  was  supported,  in  part, 
by  the  European  Research  Office  of  the  U.S.  Army,  London,  England,  under 
contract  number  DAJA37-81-C-0236. 

The  report  summarizes  briefly  the  work  detailed  in  the  previous  first 
and  second  annual  reports,  dated  May  1982  and  May  1983, respectively .  More 
emphasis  and  full  accounts  are  given  on  new  results  and  findings  in  the 
two  main  research  areas  supported  by  the  contract,  namely:  (a)  material 
characterization  of  copper  and  mild  steel  under  sequential  reversed 
torsional  impact  and  quasi-static  loading  at  ambient  temperature ,  and  (b) 
ratchetting  of  pressurized  tubes  subjected  to  radial  cyclic  thermal  gradient. 
These  results  are  original  and  may  be  regarded  as  a  new  contribution  to  the 
field  of  cyclic  behavior  of  materials  and  structures. 

The  author  would  like  to  express  his  appreciation  for  the  fruitful 
discussions  held  with  Dr.  M.M.  Megahed,  who  left  the  Department  on  a  leave 
of  absence  after  his  involvement  for  the  first  two  years  as  co-principal 
investigator  in  the  project. 

The  author  is  also  grateful  to  Mr.  M.M.  Mostafa  El-Kady,  Mr.  A.F. 
Bastawros  and  Mr.  M.A.  El-Kabbany,  graduate  students  in  the  Department,  for 
their  assistance  in  the  various  stages  of  the  experimental  work,  the  numerical 
computations  and  the  preparation  of  this  report. 
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GENERAL  INTRODUCTION 


Rigorous  understanding  of  the  effects  of  impact  loading  jn  the 
mechanical  response  of  materials  and  structures  is  essential  for  the 
optimum  design  and  safe  operation  of  many  sophisticated  engineering 
systems  and  components,  such  as  industrial  high-energy- rate  fabrica¬ 
tion  processes  and  nuclear  reactor  containments. 

Extensive  data  are  available  at  present  on  the  dynamic  behavior 
of  most  metals  in  uniaxial  tension,  compression,  torsion  and  pure  shear, 
when  they  are  subjected  to  a  variety  of  loading  conditions,  ranging  from 
those  characterised  by  monotonic  constant  strain  rates,  to  those  involv¬ 
ing  forward  or  reverse  strain- rate  jumps  of  several  orders  of  magnitude. 
What  is  missing  in  the  present  data  banks,  however,  is  detailed  informa¬ 
tion  concerning  the  mechanical  response  under  reversed  sequential  cyclic 
loading  at  impact  strain  rates.  Such  data  are  badly  needed  for  engineer¬ 
ing  design  purposes,  and  also  to  help  in  the  formulation  of  proper  con¬ 
stitutive  equations  and  the  accurate  modeling  of  deformation  processes. 

On  the  other  hand,  much  effort  has  been  directed  in  recent  years 
towards  an  understanding  of  the  mechanism  of  cyclic  inelastic  deformation 
of  structures.  Whilst  it  has  been  realized  that  purely  elastic  design  cf 
structures  is  both  wasteful  of  material  and  often  impractical,  the  presence 
of  cyclic  inelastic  strains  (ratchetting)  may  lead  to  the  failure  of 
structures  by  excessive  distortion  or  low  cycle  fatigue.  This  is  partic¬ 
ularly  the  case  in  components  subjected  to  continuous  mechanical  loading 
together  with  cyclic  thermal  loading,  such  as  the  tubular  fuel  cans  used 
in  liquid  metal  fast  breeder  reactors,  aeroplane  wings  and  off-shore 
structures.  Present  design  codes  require  that  a  design  at  specified  opera¬ 
ting  conditions  shall  be  free  from  ratchetting,  or  if  ratchetting  does 
occur,  that  the  strain  accumulation  shall  be  less  than  a  specified  limit 
during  the  expected  operating  lifetime  of  the  component. 

It  is  against  this  background  that  the  research  program  carried  out 
under  the  present  contract  was  established.  Two  main  objectives  were 
identified: 


(a)  The  novel  adaptation  of  a  well-established  impact  testing  technique  to 
the  generation  of  new  data  bases  for  OFHC  copper  and  mild  steel  regard¬ 
ing  their  behavior  under  sequential  reverse  loading.  This  is  followed 
by  an  examination  of  the  results  vis-a-vis  the  classical  plasticity 
theories,  and  the  more  realistic  constitutive  equations  recently  proposed 
in  the  literature. 

(b)  The-study  of  ratchetting  in  the  simple  but  practical  example  of  a  pressur¬ 
ized  tube  subjected  to  radial  cyclic  thermal  gradient,  as  affected  by  the 
model  used  to  describe  its  material  cyclic  behavior. 

Various  activities  were  performed  in  the  three  years  of  the  contract  to 
achieve  these  two  objectives.  These  are  detailed  in  the  previous  annual 
reports,  as  well  as  in  the  current  final  one. 

The  first  annual  report,  dated  May  1982,  presented  a  detailed  descrip¬ 
tion  of  the  modified  split  Hopkinson  torsional  bar  designed  and  manufactured, 
partly  under  the  present  contract.  The  performance,  trouble  shooting  and 
refinements  of  this  apparatus  were  also  delineated.  Moreover,  results  collec¬ 
ted  with  this  new  testing  facility  on  an  aluminium  alloy  were  presented  and 
analysed.  The  report  also  included  the  new  generalized  plane-strain  formula¬ 
tion  proposed  for  analyzing  the  case  of  a  thick  walled- cylinder  subjected  to 
a  steady  internal  pressure  and  a  cyclic  temperature  gradient.  A  review  of 
cyclic  phenomena  in  metals  indicated  that  the  kinematic  and  isotropic  harden¬ 
ing  hypotheses  provide  extreme  limits  of  material  modeling,  and  therefore  may 
be  used  in  the  analysis  of  mechanical  components  subjected  to  cyclic  loading. 
Results  obtained  with  the  kinematic  hardening  assumption  for  304  stainless 
steel  were  also  provided. 

The  second  annual  report,  dated  May  1983,  presented  a  variety  of  results 
obtained  in  the  area  of  analysis  of  structures  subjected  to  cyclic  loading. 
First,  the  problem  the  tube  under  internal  pressure  and  cyclic  non-linear 
thermal  gradient  was  investigated  using  a  uniaxial  model.  Numerical  results 
were  also  presented  for  kinematic  and  isotropic  hardening  theories  of  behavior. 
The  previously  proposed  generalized  plane-strain  formulation  was  also  used 
together  with  kinematic  hardening  assumption,  and  results  compared  for  this 
biaxial  model  with  those  of  the  uniaxial  one,  as  well  as  with  experimental 


results.  Next,  a  study  of  kinematic  hardening  analysis  of  ratchet  strain 
in  the  "Pulley  Test"  was  presented.  This  is  used  to  simulate  the  uniaxial 
model  of  a  thin  tube  subjected  to  combined  internal  pressure  and  cyclic 
thermal  stress  across  its  wall.  Finally,  the  concept  of  the  reference 
stress  method  (RSM)  and  its  applications  in  cyclic  plasticity  analysis  were 
presented. 

The  present  final  report  is  divided  in  two  parts,  and  provides  addi¬ 
tional  new  results  collected  under  the  contract.  In  Part  I,  the  current 
status  of  the  modified  split  Hopkinson  torsional  bar  is  briefly  described. 
Results  collected  for  OFHC  copper  and  mild  steel  are  then  presented  and 
discussed.  These  are  contrasted  with  predictions  of  classical  plasticity 
theories.  Recently  proposed  constitutive  equations  are  also  examined  to 
identify  those  most  probably  fitting  the  present  experimental  data.  In 
Part  II,  the  biaxial  model  of  the  thick  walled  tube  previously  studied  is 
being  considered  together  with  isotropic  hardening  theory,  and  numerical 
results  compared  with  those  previously  obtained.  Using  the  biaxial  model 
and  the  developed  numerical  method,  the  full  solution  map  for  any  specific 
problem  can  be  numerically  constructed  for  any  given  form  of  temperature 
distribution  and  any  well  mathematically-founded  hardening  rule. 


PART  I 


Behavior  of  OFHC  Copper  and  Mild  Steel  Under  Sequential  Reversed 
Torsional  Loading  at  Ambient  Temperature 


SECTION  1 

INTRODUCTION 

Conventional-type  apparatus,  in  which  the  stress  may  vary  consid¬ 
erably  along  the  length  of  the  specimen  due  to  inertia  effects,  cannot 
be  used  in  the  impact  range  of  strain  rates.  The  most  popular  experimental 
device  which  can  meet  the  high-speed  test  requirements  is  the  spiit-Hopkinson- 
bar  (SHB)  apparatus.  Over  the  past  two  decades,  the  SHB  has  been  used  repea¬ 
tedly  for  testing  materials  at  fairly  constant  strain  rates,  and  also  adapted 
successfully  to  study  history  effects  by  performing  incremental  strain-rate 
change  tests/1-3/. 

For  cyclic  loading  at  impact  strain  rates,  the  SHB  needs  further 
modifications.  For  small  plastic  straining  in  the  tension-compression  or 
torsional  mode,  the  adaptation  of  the  SHB  is  simply  based  on  the  one- 
dimensr.onal  elastic  wave  propagation  theory/4-6/.  This  can  be  done  by 
introducing  suitable  discontinuities  at  specific  locations  along  the  testing 
bar  system  which  cause  reflections  of  proper  magnitudes  in  the  travelling 
elastic  waves,  and  consequently  unloading  and  reverse  loading  in'. the  test 
specimen. 

% 

On  the  other  hand,  subjecting  a  material  specimen  in  the  SHB  to  reversed 
sequential  loading  causing  large  plastic  strains  is  not  a  straightforward 
procedure.  In  the  tension-compression  mode,  two  major  difficulties  exist. 

The  first  is  related  to  the  specimen  itself,  where  a  particular  geometry  is 
usually  required  for  the  case  of  tension  which  drastically  differs  from  that 
required  for  compression.  The  second  difficulty  relates  to  the  two  input  and 
output  bars  which  should  be  suitably  modified  to  accomodate  the  corresponding 
specimen.  Also,  in  this  t-c  mode,  only  a  single  sequence  is  possible.  Thus, 
a  specimen  can  be  first  loaded  dynamically  in  tension.  Thereafter,  compres¬ 
sion  tests  are  conducted  on  suitable  wafers  machined  from  this  previously 
strained  tensile  specimen.  This  procedure  is  accompanied  by  two  inherent 
limitations,  namely:  (a)  the  long  elapsed  time  necessary  between  loadings, 
and  which  is  usually  of  the  order  of  a  few  hours  at  least,  and  (b)  the 
relatively  small  values  of  plastic  strain  involved  to  ensure  uniform  deforma¬ 
tion,  away  from  necking  in  tension  during  prestraining,  and  from  barrelling 
in  compression  during  the  reverse  straining. 


It  seems  timely  at  this  point  to  list  some  of  the  essential  general 
requirements  needed  for  the  generation  of  relevant  data  on  material  behavior 
under  cyclic  loading  at  impact  strain  rates.  These  are:  (a)  the  development 
of  a  testing  system  capable  of  both  forward  and  reverse  loading  at  fairly 
constant  strain  rates  and  up  to  large  values  of  strain,  and  (b)  the 
capability  of  subjecting  the  same  specimen  to  load  reversals  of  small  or 
large  strain  amplitudes  within  a  reasonable  amount  of  elapsed  time.  These 
considerations  automatically  favour  the  torsion  of  thin-walled  tubes  as  the 
ideal  mode  of  loading.  This,  in  addition,  offers  the  following  advantages: 

1.  The  achievement  of  a  homogeneous  stress  state  and  a  large  range  of  uniform 
deformation  without  significant  change  in  the  specimen  shape. 

2.  The  performance  of  forward  and  reverse  loading  on  the  same  specimen,  and 
its  sequential  repetition  for  more  than  one  cycle,  up  to  fracture. 

3.  For  dynamic  testing,  the  torsional  version  of  the  SHB  has  proved  its 
superiority/2,7/,  particularly  since  many  of  the  difficulties  which  accompany 
the  axial  versions  are  eliminated. 


Satisfactory  dynamic  cycling  testing  of  materials  can  therefore  be 
accomplished  by  twisting  a  tubular  specimen  at  an  impact  strain  rate  in  a 
torsional  SHB,  up  to  a  given  strain,  and  unloading.  The  same  specimen  can 
then  be  retwisted  in  the  reverse  direction  by  simply  reversing  the  direction 
of  the  loading  torsional  wave  propagating  in  the  system  from  clockwise  to 
anti-clockwise  or  vice-versa.  This  procedure  may  then  be  repeated  again  and 
again  on  the  same  specimen  until  fracture. 

SECTION  2 

EXPERIMENTAL  DETAILS 
2. 1  Description  of  Apparatus 

In  its  final  form,  the  torsional  SHB  system  designed  and  manufactured 
under  the  present  contract  for  reversed  sequential  testing  of  materials  is 
shown  in  the  schematic  drawing  of  Fig.  1.  Design  considerations  and  complete 
description  of  each  one  of  the  system  components  as  well  as  related  instrumen 
tation,  all  in  their  preliminary  status,  have  been  presented  in  full  details 
elsewhere/8,9/.  To  avoid  repetition,  therefore,  only  a  brief  description  is 
given  here, with  emphasis  being  made  on  the  main  modifications  introduced  on 


-  8  - 


the  system  since  then,  and  on  the  main  adjustments  found  necessary  to  carry 
out  the  present  experimental  program. 

The  apparatus  essentially  consists  of  the  tvo  input  and  output  bars 
sandwiching  the  specimen,  a  frictional  clasping  device,  an  actuator 
arrangement  and  a  slow  rotational  drive  unit.  The  actuator  arrangement  is 
of  a  particularly  novel  design  and  represents  the  salient  feature  of  the 
apparatus.  It  allows  the  rotation  of  the  input  bar  in  either  the  clockwise 
or  anticlockwise  direction.  Dynamic  torsional  straining  is  achieved  by  the 
release  of  a  pre-stored  torque  pumped  by  the  actuator  into  a  portion  of  the 
input  bar.  One-half  of  this  stored  torque  travels  down  the  input  bar  to  load 
the  specimen,  while  the  other  half,  after  reflection  at' the  fixed  end  of  the 
input  bar,  unloads  it.  The  loading  incicent  wave,  characterized  by  a  short  rise 
time  and  a  conrl-iz^  amplitude  thereafter,  is  partly  reflected  and  partly  trans¬ 
mitted  at  the  specimen.  The  axial  position  of  the  frictional  clamp  controls 
the  useful  loading  duration  of  the  incident  wave,  and  consequently  determines 
the  amount  of  straining  a  specimen  may  be  subjected  to  before  being  unloaded. 

The  stored  torque  in  the  input  bar  was  measured  by  an  electrical  resis¬ 
tance  strain-gage  bridge  cemented  to  its  surface.  The  bridge  output  was  fed 
to  a  digital  hybrid  multimeter  calibrated  to  measure  torque.  On  the  other 
hand,  the  Incident  and  transmitted  waves  were  measured  by  similar  torque 
strain-gage  bridges  on  the  input  and  output  bars,  respectively.  Outputs  of 
these  two  bridges  were  fed  to  a  Tektronix  type  5111A  storage  oscilloscope, 
equipped  with  type  SA26  dual  differential  amplifier  and  type  5B12N  dual  time 
base.  Displays  on  the  oscilloscope  were  recorded  photographically  on  type 
665  Polaroid  film  using  a  type  C-5C  Tektronix  camera  (Fig.  2). 

Quasi-static  reversed  sequential  cyclic  loading  could  be  achieved  using 
the  same  bars.  Such  configuration  is  particularly  useful  when  performing 
tests  involving  a  change  in  the  strain  rate.  In  the  present  work,  however, 
a  separate  pair  of  shorter  bars  was  used.  As  shown  in  Figs.  3  and  A,  these 
bars  were  mounted  on  separate  supports  at  the  distal  end  of  the  SHB  base. 

The  inner  end  of  each' bar  was  provided  with  three  jaws  to  engage  with  similar 
ones  on  a  specimen  holder  in  order  to  facilitate  the  mounting  of  the  specimen 
assembly.  The  far  end  of  the  input  bar  was  keyed  to  a  quasi-static  rotational 
drive  unit,  while  the  far  end  of  the  output  bar  was  rigidly  fixed  in  a 


frictional  clamp.  The  drive  unit  consisted  of  a  1-HP  single-phase  electric 
motor  having  a  fixed  speed  of  1440  rpm,  two  compact  speed-reduction  units  of 
1:30  fixed  reduction  ratio  each,  a  4-speed  gearbox  of  1:2.66,  1:2,  1:1.6  and 
1:1  reduction  ratios,  and  an  electro-magnetic  brake.  These  components  are 
shown  separately  in  Fig.  5. 

During  a  test,  the  torque  and  angular  rotation  were  continuously 
plotted'on  an  x-y  recorder.  This  was  done  by  feeding  the  output  signals  of 
a  torque  strain-gage  bridge  and  a  differential  transformer  (DCDT)  to  the 
vertical  and  horizontal  channels  of  the  recorder,  respectively.  By  proper 
scaling,  this  plot  was  calibrated  to  provide  a  stress-strain  curve  for  the 
specimen  (Fig.  6).  Elastic  deformation  of  the  specimen  assembly  and  the 
loading  system  was  also  determined  and  deducted  from  the  overall  deformation 
during  each  test  in  the  strain  calculation  process. 

2.2  Materials  and  Specimens 

As  reported  previously/8,9/,  the  present  apparatus  has  been  used 
successfully  to  investigate  the  dynamic  hysteresis  loop  and  the  dynamic 
Bauschinger  effect  in  an  aluminum  alloy  corresponding  very  closely  to  American 
designation  A1  6061. 

,  < 

In  the  present  study,  two  other  materials  of  practical  importance  and 

of  different  crystallographic  structure  were  chosen, namely:  Oxygen-free  high- 

conductivity  (OFHC)  copper,  and  mild  steel.  Compositions  of  these  materials 

are  given  in  Table  I.  AFter  machining  to  the  nominal  configuration  shjwn  in 

Fig.  7,  each  specimen  was  measured  and  its  actual  dimensions  recorded.  Mild 

steel  specimens  were  tested  without  further  heat  treatment,  while  all  copper 
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specimens  were  annealed  in  a  vacuum  of  10  to  10  torr  at  370  C  for  1.5  hr. 

Temperature  within  the  furnace  did  not  vary  by  no  more  than  ±  3°C  from  one 

location  to  another.  The  furnace  was  then  turned  off  and  the  specimens 

allowed  to  coal  with  the  door  closed.  Average  grain  densities,  as  measured 

_2 

from  microscopic  examination,  were  found  to  be  161  mm  in  annealed  copper 
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specimens  and  266  mm  in  steel  specimens. 


2.3  Testing  Program 

A  main  objective  of  the  present  research  program  was  to  collect  experimen 
tal  data  on  the  sequential  reversed  loading  behavior  of  OFHC  copper  and  mild 


steel,  as  affected  by  strain  rate.  Therefore,  in  designing  the  correspond¬ 
ing  testing  progran,  the  following  facts  had  to  be  accounted  for  : 

(a)  Three  different  strain  rates  should  be  investigated,  one  quasi-static 
and  two  dynamic. 

(b)  In  using  the  SHB  apparatus  for  dynamic  torsional  testing,  it  can  be 
shown  that  the  following  simple  expressions  are  reliable  for  calculating 
the  specimen  stress,  strain  rate  and  strain  /10,7/: 


where 
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Since  the  relative  magnitudes  of  4^  and  depend  essentially  on 
the  physical  properties  of  the  material  specimen,  it  is  therefore 
expected  that,  for  a  given  specimen  geometry  and  a  given  (i.e.  fixed 
value  of  stored  torque  in  the  input  bar)  ,  upon  sequential  loading  in  the 
same  or  reverse  direction; 

(i)  the  same  strain  rate  and  total  strain  will  be  repeatedly  achieved  in 
a  non  straln*diardening  material  specimen. 

(ii)  slightly  decreasing  values  of  strain  rate  and  total  strain  amplitude 
will  sequentially  take  place  in  a  strai.  -hardening  material.  In  the 
latter  case,  however,  a  slight  incremental  increase  in  the  stored  torque 
(and  hence  in  $^)  will  keep  the  strain  rate  and  total  strain  constant 
upon  repeated  loading.  Hence,  with  the  SHB  used  in  the  present  investiga¬ 
tion,  sequential  reverse  loading  of  fixed  strain  rate  and  strain  amplitude 
can  always  be  performed. 

(c)  For  a  given  specimen  strain  rate,  it  Is  also  important  to  study  the 

effects  of  changing  the  value  of  the  fixed  strain  amplitude  characterizing 
the  series  of  reverse  sequential  impacts.  This  can  only  be  achieved  in 
the  present  apparatus  by  changing  the  total  useful  duration  of  the  loading 
pulse  through  a  physical  change  in  the  axial  location  of  the  friction 
clamp  along  the  system,  as  shown  in  Fig.  8.  Upon  reaching  the  specimen, 


the  reflected  waves  caused  by  the  fixed  hydraulic  actuator  terminate  the 
useful  duration  of  the  loading  wave  and  force  the  specimen  to  unload  from 
the  maximum  stress  already  reached. 

In  the  present  work,  three  different  locations  were  chosen  for  the  fric¬ 
tion  clamp.  These  are  designated  as  Position  I,  II  and  TII,  yielding  a 
useful  loading  time  in  the  incident  wave  of  530,  770  and  1100  ps  respec¬ 
tively.  Experimental  verifications  of  these  values  are  shown  in  the 
records*of  Figs.  9  (a),  (b)  and  (c)  corresponding  to  the  three  clamp 
positions,  respectively.  These  records  were  obtained  by  using  a  dummy 
specimen  between  the  two  bars  and  simultaneously  measuring  the  travelling 
pulses  at  the  input  bar  strain  gage  station  (top  trace)  and  the  output  bar 
strain  gage  station  (bottom  trace) . 

(d)  Changing  the  strain  rate  in  the  present  setup  can  be  achieved  simply  by 
increasing  the  amount  of  stored  torque  in  the  input  bar.  Further  coqjpli- 
oations  and  disturbances  were  found  to  accompany  this  procedure,  however, 
which  forced  the  use  of  a  constant  value,  within  the  range  of  13.3  -  14.5  kgm, 
for  the  stored  torque  in  all  dynamic  testing.  Instead,  the  strain,  where  as 
indicated  by  Eq.  (lb),  a  higher  strain  rate  is  induced  by  a  smaller  gage 
length  and  vice  versa.  Two  different  values  of  nominal  specimen  gage  length 
were  used,  namely  1.1  and  3.0  mm. 

Based  on  the  above  considerations,  the  following  characteristics  accompanied 

the  dynamic  testing  program  of  the  present  research: 

-  Stored  torque  :  13.3  -  14.5  kgm  \ 

-  Axial  position  of  cl  mp  from  actuator  end  :  75,108  and  155  cm 

-  Corresponding  useful  duration  of  incident  wave  :  530,  770  and  110  ps 

-  Nominal  specimen  gage  length  :  1.1  and  3.0  mm 

-  Type  of  test  :  Forward,  reverse,  forward,  ieverse; 

alternatively  referred  to  as  (+) ,  (-) ,  (+) ,  (-). 

-  Testing  temperature  :  22°C. 

All  specimens  used  ir.  the  present  work,  together  with  their  nominal  and 

actual  dimensions,  as  well  as  the  test  conditions  they  were  subjected  to,  are 

listed  in  Table  II. 


2.4  Testing  Procedure 

In  each  dynamic  test,  the  procedure  involved  many  steps  which  may  be 
grouped  into  three  main  phases  :  (a)  setup  preparation  involving  mainly  the 
fixation  of  the  friction  clamp  at  the  proper  axial  position,  the  mounting  of 
the  specimen  between  the  two  bars  with  Araldite  epoxy  adhesive,  and  the  adjus¬ 
tment  of  the  hydraulic  actuator;  (b)  instrumentation  checkup;  and  (c)  specimen 
loading.  The  two  last  phases  are  repeated  on  the  same  specimen  for  further 
loading  in  the  reverse  direction,  according  to  the  previous  program.  The 
whole  testing  procedure  is  summarized  in  the  flow  chart  of  Fig.  10. 

Quasi-static  testing  involved  the  separate  fixation  of  the  specimen  in 
its  holder  with  Araldite  epoxy  adhesive,  and  the  fixation  of  this  assembly 
to  the  pair  of  jaws  of  the  two  short  input  and  output  bars.  The  distal  end  of 
the  output  bar  is  then  clamped,  the  instrumentation  checked  and  the  electric 
motor  started  to  load  the  specimen  at  the  required  slow  speed  to  a  pre-selected 
plastic  strain  in  the  forward  direction.  The  clamp  is  then  opened  to  release 
the  stored  torque  in  the  system,  re-tightened,  and  the  torque  re-applied  at  the 
same  rotational  speed  to  the  same  value  of  plastic  strain  in  the  reverse  direc¬ 
tion.  This  procedure  is  then  repeated  any  number  of  times  as  required,  with 
the  torque  versus  angle  of  rotation  at  the  motor  end  of  the  input  bar  being 
plotted  on  the  x-y  recorder,  in  each  case. 

2 . 5  Typical  Test  Records 
2.5.1  Quasi-static  loading 

Figures  11  and  12  represent  typical  torque-twist  angle  records  obtained 
on  copper  and  steel,  respectively.  Each  record  was  obtained  on  a  separate 
specimen  which  was  loaded  in  the  forward  direction  up  to  a  given  value  of 
strain,  unloaded  by  the  release  of  the  friction  clamp,  then  immediately  reloaded 
at  the  same  twisting  rate  in  the  opposite  reverse  direction.  Eiples  exhibited 
on  the  traces  do  not  indicate  intrinsic  yield  point  behavior  or  material 
instabilities,  but  are  essentially  due  to  the  relatively  high  amplification 
■jsed  on  the  x-y  recorder.  Traces  were  smoothed  out  upon  data  reduction. 


Forward  yielding  of  copper  is  very  much  undefined  and  seems  to  occur  at 
ery  low  torque.  The  reverse  loading  is  characterised  by  a  much  higher  slope 
f  the  elastic  line  and  a  sharper  yield  point  than  in  forward  loading.  In 
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contrast,  the  forward  yield  of  steel  Is  more  well  defined,  while  the  elastic 
behaviour  has  almost  the  same  slope  in  reverse  and  forward  loading. 

2.5.2  Dynamic  loading 

A  representative  collection  of  dynamic  test  records,  reduced  in  scale 
for  ease  in  comparison,  is  presented  in  the  matrix  of  Pigs.  13  and  14  for 
copper  and  steel,  respectively.  Each  matrix  consists  of  24  records.  Any 
four  records  gathered  together  in  one  row  are  the  results  of  testing  on  a 
single  specimen  in  the  following  sequence:  forward,  reverse,  forward,  reverse. 
The  first  three  rows  correspond  to  dynamic  testing  at  the  same  impact  strain 
rate  y^,  while  they  differ  among  themselves  in  the  value  of  the  constant  strain 
amplitude,  yq,  achieved  in  each  of  their  four  sequential  loadings.  These  same 
three  values  of  strain  amplitude  were  also  achieved  in  the  tests  represented 
by  the  records  of  the  second  set  of  three  rows;  the  Impact  strain  rate  for  this 
set,  however,  being  smaller  than  y^ 

In  all  oscillograms,  the  upper  trace  shows  the  input  bar  gage  sigr  ul  (i.e. 
the  incident  and  reflected  tcrque),  while  the  lower  trace  shows  the  output-bar 
gage  signal  (i.e.  output  torque).  For  all  traces  in  both  Figs.  13  and  14,  the 
horizontal  sweep  was  kept  the  same  at  200  ps/div.  The  veritcal  sensitivity  for 
the  upper  traces  shown  in  the  four  records  in  each  row  were  identical,  but 
smaller  than  that  used  for  the  lower  traces  in  the  same  row.  Traces  of  corres¬ 
ponding  records  in  different  rows,  however,  are  not  directly  comparable,  since 
the  amplitudes  obtained,  besides  being  function  of  the  oscilloscope  settings, 
depend  also  onthe  excitation  voltage  of  the  strain  gage  bridges  which,  being 
provided  by  dry  batteries,  varied  from  one  test  to  the  other.  Table  III  gives 
in  detail  the  calibration  of  the  vertical  traces  in  all  records.  It  should  also 
be  noted  in  Figs.  i3  and  14  that,  before  each  impact  in  the  loading  sequence, 
the  positions  of  the  baseline  of  the  top  and  bottom  traces  were  changed  in  order 
to  accomodate  the  expected  trace  motion. 

0 

2. 6  Data  Reduction 
2.6.1  Quasi-static  response 

In  quasi-static  testing,  the  elastic  response  of  the  system,  apart  from 
the  specimen,  was  drawn  as  a  straight  line  on  the  same  torque-twist  diagram 
obtained,  whether  in  the  forward  or  reverse  loading,  in  order  to  give  the 


elastic-plastic  specimen  behavior.  Shear  stress  was  calculated  from  the  relation: 


t  -  T/2  ir  r2t  ,  (3) 

s 

whereas  shear  strains  were  calculated  by  measuring,  at  various  values  of  the 
torque,  offset  twist  values  between  the  experimental  total  torque-twist  plot  and 
using  the  equation  : 


Y*  0r/l  .  (4) 

S 

It  should  be  noted  that  throughout  this  report,  "engineering"  shear  strains  are 
quoted,  l.e.  twice  the  tensor  values.  .■ 

2.6.2  Dynamic  response 

In  this  case,  oscilloscope  records  in  each  impact  loading  were  analysed 
in  terms  of  specimen  stress,  strain  and  strain  rate  following  the  usual  technique 

•  i  *  ‘ 

for  a  SHB  test.  This  is  based  on  the  assumption  of  uniform  stress  state  in  the 
short  specimen  gage  length,  which  can  be  justified  in  particular  when  a  torsional 
pulse  of  constant  amplitude  is  applied  as  in  the  present  apparatus . 

Specimen  stress-time  behavior  is  usually  determined  from  the  transmitted 
pulse  as  measured  on  the  output  bar,  while  the  strain  rate-time  characteristic, 
and  consequently  strain  by  Integration,  can  be  obtained  from  the  reflected  pulse 
if  this  can  be  completely  measured  on  the  input  bar.  As  shown  in  Fig.  1,  the 
strain-gage  bridge  or.  the  input  bar  is  positioned  far  enough  from  the  specimen, 
so  that  a  large  proportion  of  the  incident  pulse  can  be  measured  and  judged  in 
quality  before  the  reflected  pulse  from  the  specimen  arrives  at  the  gages.  This 
Incident  pulse  was  always  found  to  be  of  constant  amplitude  immediately  following 
a  short  rise  time,  as  can  be  seen  from  the  oscilloscope  traces  shown  previously 
in  Figs.  9(a)-(c).  Therefore,  the  amplitudes  of  the  reflected  pulse  were  calcula¬ 
ted  as  the  difference  between  the  absolute  values  of  the  incident  and  transmitted 
pulses  taken  at  coincident  times.  This  eliminated  the  necessity  of  a  separate 
direct  measurement  of  the  reflected  pulse  with  another  set  of  strain  gages. 

The  analysis  of  the  test  records  was  carried  out  by  using  the  standard 
expressions  given. by. Eq*._(l) ,  relating.. shear. stress ,  shear. strain  rate  and 
shear  strain  to  the  Incident  and  transmitted  pulses;  calibration  being  effected 
by  using  shunt  resistances  across  the  strain  gauges/7/. 


All  calculations  were  performed  on  a  Sinclair  ZX  Spectrum  desk-top 
calculator  equipped  with  an  Alpha  32  printer.  The  computer  program  used  is 
very  similar  to  that  used  previously  and  listed  elsewhere/ 11/. 

SECTION  3 

RESULTS  AND  DISCUSSIONS 

3. 1  Copper 

3.1.1  Sequential  reverse  loading  response 

The  stress-strain  characteristics  of  OFHC  copper  specimens  of  1.1  mm 
nominal  gage  length  under  sequential  reverse  twisting  at  the  strain  rate  of 
about  1100  s  *  are  presented  in  Figs.  15(a),  (b)  and  (c)  for  the  cyclic 
strain  amplitude  of  about  0.46,  0.68  and  0.97,  corresponding  to  clamp  position  I, 

II  and  III,  respectively.  Four  loading  sequences  are  shown  in  each  of  the  three 
cases,  causing  a  total  cumulative  strain  of  about  1.32,  1.75  and  3.7,  respectively, 
after  which  the  test  was  stopped  although  the  specimen  did  not  show  any  sign  of 
failure  of  flow  instability.  The  elapsed  time  between  any  two  consecutive  load¬ 
ings  did  not  exceed  five  minutes  in  all  cases.  Note  that  the  curves  resulting 
from  reverse  loadings  (-)  are  plotted  in  the  figures  upward,  instead  of  downward, 
for  ease  of  comparison  with  the  curves  resulting  from  forward  loadings  (+) .  As 
discussed  above,  it  is  difficult  with  strain-hardening  materials  to  achieve  the 
same  strain  rate  by  using  the  same  amount  of  stored  torque  in  consecutive  load¬ 
ings.  A  decreasing  strain  rate  is  obtained  instead,  as  observed  in  Figs.  15(a) 
and  (b) .  A  slight  increase  in  the  amount  of  stored  torque,  however,  can  help  in 
keeping  the  strain  rate  constant,  as  seen  in  Fig.  15(c).  The  effects  of  this 
decrease  in  strain  rate  between  consecutive  loadings  were  judged  to  be  of ^secondary 
importance  and,  therefore,  neglected  in  the  present  analysis. 

Close  examination  of  the  response  reveals  the  following  behavioral  features: 
(a)  The  first  loading  is  characterized  by  an  elastic  behavior  of  small  slope  and 
an  undefined  transition  to  plastic  flow.  Yield  stress  level  was  taken  at  the 


limit  of  proportionality,  and  is  pointed  at  in  the  figures  by  a  short  horizon¬ 
tal  arrow.  Following  this  initial  yield,  plastic  deformation  takes  place 
with  considerable  work  hardening,  as  it  is  well  known  for  annealed  pure  copper. 
This  is  described  by  an  equation  of  the  form: 


where  K  and  n  are  the  strain-hardening  coefficient  and  exponent, 
respectively,  whose  calculated  values  are  shown  in  Table  IV.  Eq.  ^5) 
was  assumed  to  be  valid  for  large  strains,  and  was  consequently  used 
to  extrapolate  the  flow  curve  beyond  the  level  reached  upon  first  load¬ 
ing,  as  shown  by  dotted  lines  in  Figs.  15(a),  (b)  and  (c). 

(b)  Sequential  reverse  loadings  are  characterized  by  initial  elastic 
behavior  with  a  slope  considerably  larger  than  that  of  first  loading 

and  unloading.  The  limit  of  proportionality  is  taken  as  the  yield  point, 
whose  level  is  always  found  lower  than  the  maximum  stress  level  reached 
in  the  previous  loading,  indicative  of  the  existence  of  a  Bauschinger 
effect. 

(c)  In  Fig.  15(a),  where  the  strain  amplitude  is  the  smallest,  the  flow  curve 
upon  sequential  reverse  loading  starts  by  being  higher  than  the  level 
expected  by  extrapolation  from  continuous  forward  loading,  but  subsequently 
gets  closer  to  this  extrapolated  level,  and  ultimately  ends  up  below  it  in 
the  fourth  loading,  indicative  of  continuous  "softening"  of  the  material 
with  reverse  loading.  Moreover,  the  flow  curve  flattens  out  with  increasing 
sequential  loading,  probably  reaching  saturation  at  a  later  stage.  This 
same  trend  is  also  seen  to  occur,  but  more  rapidly,  in  the  case  of  the  larger 
strain  amplitude  of  Fig.  15(b),  whereas  "softening"  starts  from  the  first 
reverse  loading  in  the  case  of  the  largest  strain  amplitude  of  Fig.  15(c). 

In  all  three  cases,  it  also  appears  that,  irrespective  of  the  previous  load¬ 
ing  history,  the  flow  curve  in  sequential  reverse  loading  assumes  a  lower 
level  than  that  of  continuous  forward  loading  starting  from  a  total  cumulative 
strain  of  about  0.95. 

The  stress-strain  characteristics  of  OFHC  copper  specimens  of  3.0  mm 
nominal  gage  length  under  sequential  reverse  twisting  at  the  strain  fate  of 
about  550  s  are  presented  in  Figs.  16(a),  (b)  and  (c)  for  the  cyclic  strain 
amplitude  of  0.23,  0.33  and  0.46,  corresponding  to  clamp  position  I,  II  and  III, 
respectively.  Exactly  the  same  characteristics  and  behavioral  features  are 
exhibited  here,  as  those  discussed  above  pertaining  to  the  higher  strain  rate  of 
1100  s  1.  In  the  present  case,  however,  the  flow  curve  in  sequential  reverse 
loading  assumes  a  lower  level  than  that  of  continuous  forward  loading  starting 
from  a  total  cumulative  strain  of  about  0.60. 


’>  -JI  ■-» 
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At  the  quasi-static  strain  rate  of  about  0.001  s  ,  the  stress-strain 
characteristics  of  OFHC  upon  forward  and  reverse  twisting  for  the  prestrain 
of  0.43,  0.76  and  1.1  are  presented  in  Figs.  17(a),  (b)  and  (c)  respectively. 
Similar  characteristics  are  exhibited  with  respect  to  forward  loading.  Upon 
reverse  loading  at  the  smallest  prestrain,  the  flow  curve  also  assumes  a 
higher  level  than  that  of  the  continuous  forward  loading  obtained  by  extra¬ 
polation,  but  the  difference  is  not  large.  For  larger  values  of  prestraining, 
the  reverse  flow  curve  is  characterized  after  the  transition  to  plastic  flow 
by  an  almost  constant  stress  level  plateau,  which  lasts  for  longer  reverse 
deformations,  the  larger  the  forward  prestrain.  This  plateau  is  then  followed 
by  strain  hardening,  but  the  flow  curve  level  remains  lower  than  that  of  con¬ 
tinuous  forward  loading. 

At  all  strain  rates,  therefore,  sequential  reverse  loading  causes  the 
flow  curve  to  flatten  out  and  to  ultimately  reach  saturation.  This  can  be 
attributed,  as  suggested  by  Armstrong  et  al/12/  who  exposed  annealed  1100 
aluminum  specimens  to  slow  multiaxial  straining,  to  the  large  accumulated  plastic 
strain  and  to  an  enhancement  of  ’softening’  or  "recovery"  processes. 

3.1.2  Analysis  of  behavior  in  terms  of  B.E.  parameters 

Some  of  the  features  discussed  above  may  now  be  described  quantitatively, 

as  shown  in  Fig.  IE  defining  different  parameters  which  may  be  used  for  that 

purpose.  At  *  is  the  amount  of  strain  hardening  during  prestraining  . 

t  is  the  reverse  yield  stress,  and  x_  «  t  -  It  I  is  the  decrease  of  yield 
yr  Bo*  yr 1 

stress,  commonly  denoted  as  "Bauschinger  stress".  In  the  present  analysis, 
the  magnitude  of  the  Bauschinger  effect  is  expressed  by  two  indicators:  T 

B 

and  T  /At,  the  latter  parameter  being  more  important  from  the  viewpoint  of 
microscopic  aspects  of  the  phenomenon/ 13/.  On  this  basis^  Table  V  lists  various 
quantities  as  calculated  from  the  first  reversal  in  Figs.  15,  16  and  17. 

The  reverse  yield  stress,  t  ,  and  the  Bauschinger  stress,  t'_,  are  both 

yr  B 

plotted  versus  the  prestrain  y  ,  for  various  strain  rates,  in  Fig.  19.  By 
definition,  x  should  equal  T  and  take  the  values  40,  25  and  20  MPa  for  the 

yr  y  -i 

three  strain  rates  1100,  550  and  G.001  s  ,  respectively, at  Y  »  Y  -  0,  while 

o  y 

x„(  *  t  -  t _ )  should  always  vanish  at  y  -  0;  such  trends  are  exhibited  to 

d  o  yr  o 

some  extent  if  the  shown  curves  are  extrapolated,  as  shown  by  the  dashed  lines, 
to  zero  prestrain.  As  can  be  seen  from  Fig.  19,  for  a  given  y»  increasing  YQ 


results  in  an  Increase  in  t  and  in  t_.  On  the  other  hand,  at  any  given  y  > 

yr  o  o 

increasing  y  increases  r  and 

yr  B 

Figure  20  shows  the  effect  of  prestrain  and  strain  rate  on  the  second 
Bauschinger  effect  Indicator  t./At.  This  ratio  Increases  with  increasing 
prestrain,  and  decreases  with  increasing  the  impact  strain  rate.  The  present 
results  are  obtained  at  large  values  of  prestrain,  in  contrast  to  those  avail¬ 
able  in  the  literature  for  copper  which  are  at  much  lower  prestrains  / 12/ ; 
moreover, *the  definition  used  for  the  yield  stress  in  materials  like  copper 
seems  to  affect  the  trend  of  behavior.  On  the  other  hand,  comparing  the 
present  results  with  those  obtained  previously  for  an  aluminum  alloy  /9/  does 
not  indicate  any  close  correspondence,  although  the  two  materials  have  an  FCC 
crystallographic  structure.  This  essentially  indicates  the  effect  of  other 
structural  parameters,  particularly  the  stacking-fault  energy  (SFE),  on  the 

behavior  of  metals,  as  discussed  by  previous  investigators  /12,13,14/; 

•  »  * 

aluminum  being  known  to  have  high  SFE,  whereas  copper  is  known  to  have  a  low  SFE. 

The  effect  of  sequential  reversal  of  loading  is  also  examined  by  calcula¬ 
ting  the  Bauschinger  effect  parameters  corresponding  to  the  fourth  loading,  and 
contrasting  them  with  those  presented  and  studied  above  for  the  second  loading. 
All  relevant  quantities  are  listed  in  Table  VI,  as  calculated  from  the  data 
shown  in  Figs.  15  and  16.  The  reverse  yield  stress  is  plotted  versus  the 
cumulative  strain  yc»  for  the  two  impact  strain  rates,  in  Fig.  21.  Shown  in 
solid  lines  are  the  values  calculated  for  the  second  loading,  while  in  dashed 
lines  are  the  values  corresponding  to  the  fourth  loading.  For  each  strain  rate, 
values  of  seem  to  fall  on  the  same  continuous  curve  which  tends  to  a  steady- 
state  or  saturation  value  with  increasing  cumulative  shear  strain;  this  satura¬ 
tion  value,  is  higher  the  higher  the  strain  rate,  and  is  also  expected  to  be  a 
function  of  the  value  of  the  cyclic  strain  amplitude.  Such  trend  falls  within 
the  general  observations  associated  with  fatigue  and  multidirectional  testing  of 
metals  / 12,15—17/  at  low  strain  rates. 

Also  in  Fig.  21,  the  variation  of  the  Bauschinger  stress  at  the  three 

strain  rates  with  cumulative  strain,  y  ,  is  shown  in  solid  lines  for  the  second 

c 

loading,  and  in  dashed  lines  for  the  fourth  loading.  It  should  be  noted  here 

that,  in  both  cases,  *  t  -  t  where  t  is  the  maximum  stress  level  reached 

.  _  a  o  yr  o 

on  the  immediate  previous  loading.  With  both  t  and  t  tending  to  almost  steady 
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state  values  at  large  cumulative  strains,  the  same  trend  Is  also  approximately 

observed  for  T_.  Similar  observations  are  also  made  from  Fig.  22  for  the 
B 

variation  of  the  B.E.  indicator,  Tg/Ax,  with  where  At  *  tq  -  Ty. 

3.2  Mild  Steel 

3.2.1  Sequential  reverse  loading  response 

The  stress-strain  characteristics  of  mild  steel  specimens  of  1.1  mm 
nominal  gage  length  under  sequential  reverse  twisting  at  the  strain  rate  of 
about  760  s”^  are  presented  in  Figs.  23(a),  (b)  and  (c)  for  the  cyclic  ampli¬ 
tude  of  about  0.262,  0.46  and  0.65,  corresponding  to  clamp  position  I,  II  and 
III,  respectively.  Four  loading  sequences  are  shown  in  each  case,  except  in 
the  case  of  cyclic  strain  amplitude  of  0.262  shown  in  Fig.  23(a)  where  the 
test  was  ended  only  after  two  loadings  due  to  the  fracture  of  the  specimen 
flange- to-input  bar  epoxy  joint.  Total  cumulative  strains  of  about  i.7  and 

2.1  were  achieved  in  the  last  two  cases,  respectively,  and  caused  the  specimen 
to  fracture.  The  elapsed  time  between  any  two  consecutive  loadings,  like  with 
copper,  did  not  exceed  five  minutes  in  all  cases.  Note  that  here  also,  the 
curves  resulting  from  reverse  loadings  (-)  are  plotted  in  the  figures  upward, 
instead  of  downward,  for  ease  of  comparison  with  the  curves  resulting  from 
forward  loadings  (+).  Because  strain  hardening  is  quite  "mall  in  mild  steel 
in  comparison  to  copper,  the  strain  rate  can  be  seen  to  vary  very  little  upon 
sequential  loading. 

Close  examination  of  the  response  reveals  the  following  behavioral 
features: 

(a)  The  first  loading  is  characterized  by  an  elastic  behavior  of  a  large 

slope  as  expected.  Yield  stress  was  taken  at  the  limit  of  proportion¬ 
ality,  and  is  pointed  at  in  the  figures  by  a  short  horizontal  arrow. 
Following  yielding,  plastic  deformation  takes  place  and  work  hardening 
quickly  assumes  its  steady-state  value,  which  is  quite  smaller  than 
for  copper.  A  power  law  can  also  be  used  to  describe  plastic  flow  in 
mild  steel.  Corresponding  calculated  values  of  K  and  n  in  Eq.  (5)  are 
listed  in  Table  VII.  This  equation  was  assumed  to  be  valid  for  large 
strains,  and  was  used  to  extrapolate  the  flow  curve  beyond  the  level 
reached  upon  first  loading,  as  shown  by  dotted  lines  in  Figs.  23(a), 

(b)  and  (c). 


(b)  Sequential  reverse  loadings  are  characterized  by  initial  elastic 
behavior  with  a  slope  almost  similar  to  that  in  first  loading  and 
unloading.  The  limit  of  proportionality,  also  pointed  at  with  a  short 
arrow,  was  taken  as  the  yield  point. 

(c)  Plastic  flow  following  reverse  yielding  proceeds  always  at  a  level  below 
that  expected  by  extrapolation  from  continuous  forward  loading.  This 
level  gets  lower  and  lower  with  subsequent  reverse  loading,  indicative  of 
continuous  softening  of  the  material.  Also,  the  flow  curve  flattens  out 
with  increasing  sequential  loading,  most  probably  reaching  saturation  had 
the  specimen  not  fractured. 

The  stress-strain  characteristics  of  mild  steel  specimens  of  3.0  mm 
nominal  gage  length  under  sequential  reverse  twisting  at  the  strain  rate  of 
350  s  ^  are  presented  in  Figs.  24(a),  (b)  and  (c)  for  the  cyclic  strain 
amplitude  of  0.15,  0.25  and  0.312,  corresponding  to  clamp  position  I,  II  and 
III,  respectively.  Exactly  the  same  characteristics  are  exhibited  here,  as 
those  discussed  above  pertaining  to  the  higher  strain  rate. 

At  the  quasi-static  strain  rate  of  about  0.001  s  \  the  stress-strain 
characteristics  upon  forward  and  reverse  twisting  for  the  prestrains  of  0.135 
and  0.26  are  presented  in  Figs.  25(a)  and  (b),  respectively.  Similar  features 
as  those  discussed  above  for  the  impact  strain  rates  are  also  exhibited  in  the 
present  case. 

At  all  strain  rates,  therefore,  sequential  reverse  loading  causes  the 
flow  curve  to  continuously  deviate  from  the  unidirectional  continuous  flow 
curve,  and  to  flatten  out  progressively  reaching  saturation. 

3.2.2  Analysis  of  behavior  in  terms  of  B.E.  parameters 


Table  VIII  lists  various  quantities  describing  the  Bauschinger  effect 
phenomenon,  as  calculated  from  the  first  reversal  in  Figs.  23,  24  and  25. 

The  reverse  yield  stress  and  the  Bauschinger  stress  are  both  plotted  versus 
prestrain  for  the  three  strain  rates  in  Fig.  26.  Both  quantities  increase  with 
increasing  strain  rate,  at  any  given  prestrain.  Thus,  the  loss  of  strength  is 
greater,  the  higher  the  strain  rate.  This  is  also  greater  the  larger  the  pre¬ 
strain.  The  same  trends  are  also  substantiated  in  the  plots  of  Fig.  27  showing 
the  variation  of  the  second  'B.E.  parameter  t«/At  with  prestrain. 


To  contrast  these  results  with  those  corresponding  to  the  fourth  sequen¬ 
tial  loading,  all  relevant  quantities  were  calculated  for  both  cases  from 
•Figs.  23  and  24  and  listed  in  Table  IX.  and  x^  are  plotted  in  Fig.  28, 
while  x_/Ax  is  plotted  in  Fig.  29,  versus  cumulative  shear  strain,  for  the 

D 

two  impact  strain  rates,  as  solid  lines  for  the  second  loading  and  as  dashed 

lines  for  the  fourth  loading.  A  steady-state  or  saturation  value  seems  to  exist 

for  T  ,  of  a  level  being  higher  for  the  larger  strain  rate.  This  is  also 

dependent  on  the  value  of  the  cyclic  strain  amplitude  but  to  a  lesser  extent 

than  with  copper  due  to  the  difference  in  SFE  and  consequently  in  the  recovery 

processes  involved,  as  has  been  suggested  in  the  literature  by  comparing  the 

behavior  of  copper  and  iron  under  fatigue  testing  at  low  strain  rates  / 12/ . 

x_  and  x  /Ax,  used  in  the  present  context  as  Bauschinger  effect  indicators, 

B  B 

exhibit  a  similar  trend  of  behavior  as  function  of  cumulative  strain.  Both 
increase  with  increasing  sequential  loading,  but  seem  to  stabilize  at  large 
values  of  strain  .  Comparing  also  the  behavior  exhibited  in  Fig.  29  with  that 
for  copper  in  Fig.  22,  it  is  clear  from  the  shown  values  of  x  /Ax  at  a  given 

D 

Y  and  the  lowest  impact  strain  rate  that  the  Bauschinger  effect  phenomenon  is 
c 

more  appreciable  in  mild  steel  than  in  copper. 

From  the  present  results  of  both  copper  and  mild  steel,  as  well  as  from 
those  previously  reported  for  aluminum  /9/,  it  is  possible  to  conclude  that 
sequential  reverse  loading  at  large  plastic  strains  progressively  weakens  the 
material,  more  seriously  so  at  impact  strain  rates  than  at  the  quasi-static 
ones.  These  results  should  be  accounted  for  in  the  design  of  structural  engineer 
ing  components,  as  wellas  in  the  planning  of  metal  forming  operations. 

3.3  Application  of  Classical  Anisotropic  Hardening  Theories  to  Results 

It  is  clear  from  the  results  presented  above  that  the  flow  stress  of 
copper  and  mild  steel  at  quasi-static  and  at  impact  rates  of  straining  is 
direction  dependent.  This  fact  has  been  known  for  long  for  many  materials 
subjected  to  a  non-proportionate  loading  path,  at  low  strain  rates.  Various 
theories  have  been  proposed  to  account  for  this  deformational  anisotropy.  Those 
used  in  the  present  comparison  are  listed  and  discussed  briefly  in  the  following. 
Bauschinger  effect  in  shear  is  also  defined  as  predicted  in  each  case  for  later 
comparison  with  the  present  experimental  results.  Further  details  of  derivations 
can  be  found  elsewhere  /18,19/. 
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3.3.1  Prager’s  kinematic  hardening  rule  / 20/ 

Here,  linear  work  hardening  is  assumed  with  the  initial  yield  surface 

translating  rigidly  in  the  direction  of  the  exterior  normal  at  the  stress 

point.  This  predicts  that  the  difference  between  and  t  which  defines 

the  elastic  range,  upon  unloading  and  reverse  loading  from  r  ,  is  a  constant 

o 

2  x  ,  cf.  Fig.  9. 

y 

3.3.2  Ziegler* s  rule  /2 1/ 

The  proposal  here  is  that  the  initial  yield  surface  translates  rigidly 
in  the  direction  of  the  vector  connecting  its  centre  to  the  stress  point  on 
it.  This,  together  with  the  linear  work  hardening  law,  predict  an  identical 
Bauschinger  effect  in  torsion  as  in  Prager's  rule. 

3.3.3  Baltov  and  Sawczuk  / 22/ 

There  authors  used  the  kinematic  hardening  rule  to  develop  a  hardening 
model  which  allows  an  initial  von  Mises  yield  surface  to  translate,  expand 
and  rotate  in  stress  space.  A  hardening  constant  H  indicates  whether  the 
surface  expands  (H  <  0)  or  contracts  (H  >  0).  With  linear  work-hardening  , 
and  proper  boundary  conditions,  the  Bauschinger  effect  in  torsion  associated 
with  this  rule  is  given  by 

/ (-^  -  l)2  +  H  -  l)4  -  1 

ize.Zo  r  y  y  _  |  . 

Ty  Ty  /  (1°  _  d2  +  H  (-  -  i)4  +  i 
Ty  Ty 

When  H  =  0,  Eq.  (6)  reduces  to  kinematic  hardening. 


3.3.4  Edelman  and  Drucker  /23 / 

Here,  an  anisotropic  yield  function  is  developed  which  is  based  upon 
the  interaction  of  prestrain  and  yield  stress  components.  The  associated 
Bauschinger  effect  is  given  by 


T 


yr 


), 


(7) 


where  m  is  an  empirical  constant  describing  the  translation  of  the  initial 

p 

yield  locus,  and  is  the  plastic  prestrain. 
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3.3.5  Comparison  with  experimental  results 

Although  the  present  data  are  collected  mostly  at  high  strains  and  strain 
rates,  which  are  not  been  accounted  for  in  the  classical  plasticity  theories 
where  deformations  are  assumed  to  be  small,  slow  and  isothermal,  it  seems 
worthwhile  to  make  the  comparison  and  to  appreciate  the  deviations,  if  any. 
Consequently,  the  values  of  x ■  and  x^  have  been  normalized  with  respect  to  x^ 

and  listed,  together  with  previous  data  in  Tables  V  and  VIII  for  copper  and 
mild  steel,  respectively. 

Values  of  t  /x  are  plotted  versus  x  /x  in  Fig.  30  for  the  present  data 

yr  y  o  y 

on  copper  and  mild  steel,  and  also  for  the  previously  reported  data  on  aluminum  /9/. 

Various  lines  are  shown  in  this  figure.  One  of  these  lines  corresponds  to 

isotropic  hardening,  for  which  x  /x^  »  -  xq/x^;  this  is  of  course  unrealistic 

with  the  data  points.  Another  line  has  the  equation:  x  /x  *  (x  /x  )  -  2, 

yr  y  o  .  y 

which  represents  Prager’s  and  Ziegler’s  kinematic  hardening  rule;  as  discussed 
above,  this  is  also  represented  by  Eq.  (6)  with  H  *  0.  Kinematic  hardening  is 
seen  to  be  an  improvement,  but  underestimates  the  magnitude  of  the  Bauschinger 
effect.  Eq.  (6)  representing  Baltov  and  Sawczuk's  rule  is  plotted  also  in 
Fig.  30  for  different  values  of  the  hardening  constant  H.  Some  further  improve¬ 
ment  seems  to  be  introduced  in  the  predictions  of  the  relation  between  Tyt/Ty 
and  xq/x^,  but  still  no  complete  agreement  seems  to  exist,  as  expected, 
particularly  for  copper  whose  work  hardening  characteristic  is  well  rounded 
and  cannot  be  expressed  by  a  linear  law.  Reasonable  agreement  is  noticed  with 
aluminum  for  a  value  of  H  -  10  -  100.  Interestingly  enough,  this  material  was 
found  to  be  almost  insensitive  to  strain  rate  and  to  have  an  approximately 
linear  strain  hardening  /9/.  On  the  other  hand,  the  trend  of  the  agreement 
between  data  points  and  theory  seems  to  indicate  that  the  values  of  H  correspond¬ 
ing  to  the  three  materials  are  positive,  which  qualitatively  indicate  the 
Bauschinger  effect  corresponds  to  a  contraction  and  a  translation  of  the  Initial 
yield  locus.  Similar  overall  conclusions  have  been  also  reported  recently  for 
quasi-static  data  on  a  mild  steel  and  two  aluminum  alloys  / 24/ • 

Equation  (7)  can  be  rewritten  in  normalized  form  as 


x  +  x 
m  _  yr  o 

xw  •  T  ■ 


v  P 
rO 


(8) 


Checking  the  validity  of  the  Bauschinger  constant  m  in  the  theory  of  Edelman 
and  Drucker  can  be  done  by  calculating  m/xy  from  the  data  listed  in  Tables  V 


and  VIII  for  copper  and  mild  steel,  and  from  previous  data  on  aluminum  /9/» 

p 

and  plotting  m/t^  versus  Yq  •  This  is  shown  for  the  data  collected  at  various 
strain  rates  for  the  three  materials  in  Fig.  31.  For  all  three  materials  and 


all  rates,  it  is  clear  that  a  large  variation  exists  in  m/ty*  indicating  that 


Bauschinger  effect  cannot  be  controlled  by  a  single  constant  m  as  in  Edelman 
and  Drucker's  theory.  Other  parameters  and  functions  may  be  used  as  implied 
by  the  modifications  introduced  by  Yoshimura  /25 /,  Williams  and  Svensson  /26/ 
as  well  as  others,  but  these  are  not  discussed  further  in  this  report.  The 
main  conclusion  is  that  theories  of  work  hardening  based  on  classical  plasticity 
concepts  are  not  capable  of  describing  faithfully  the  reverse  loading  character¬ 
istics  exhibited  by  metallic  materials,  particularly  when  rate  effects  are  also 
predominant. 


3.4  Identification  from  the  Literature  of  Constitutive  Equations  Suitable  for 
Modeling  Bauschinger  Effect  in' Metals 
As  illustrated  from  the  experimental  results,  the  mechanical  response  of 
aluminum,  copper  and  mild  steel  under  the  present  testing  conditions  of  sequen¬ 
tial  reverse  twisting  at  large  strains  and  quasi-static  as  well  as  impact  strain 
rates  has  been  found  to  be  generally  strain-history  dependent  as  well  as  sensi¬ 
tive  to  various  degrees  to  strain  rate.  The  evaluation  of  this  response  in  the 
light  of  current  constitutive  equations  for  metals  is  not  straightforward  nor 
simple,  as  will  be  soon  appreciated. 

Presently  available  constitutive  models  for  predicting  inelastic  behavior 
in  metals  are  quite  diverse  but  it  seems  possible  to  place  them  into  a  common 
framework.  Accordingly  /26 /,  these  theories  can  be  grouped  into  one  of  three 
general  categories,  namely: 

(a)  Extensions  of  classical  rate-independent  plasticity  theories  which  are  based 
on  the  concept  of  a  yield  function.  These  can  be  divided  into:  (i)  the  un¬ 
coupled  theories  in  which  the  inelastic  deformation  is  assumed  to  be 
decomposable  into  rate-dependent  creep  and  rate-independent  plasticity 
components  /27t29/,  and  (ii)  the  unified  theories,  wherein  all  inelastic 
deformation  is  embodied  in  a  single  term  /30-33/. 

(b)  Non-linear  viscoelasticity  theories  which  are  usually  based  on  macrophenomen- 
ological  thermodynamics  or  resulting  mechanical  analogs  734-38/ • 


(c)  Microphenomenological  theories  which  find  their  origin  in  microphysics  / 39-46/ . 

Of  all  these  theories,  those  falling  in  the  last  category  have  been  found 
more  feasible  since,  in  the  most  part,  they  are  fundamentally  based  in  the 

micromechanibBS  causing  inelastic  deformation,  such  as  dislocation  glide  and 
climb,  dislocation  pileups,  cross  slip.  grain-boundary  sliding,  diffusion,  etc. 
Consequently,  many  of  these  theories  were  subjected  recently  to  considerable 
study,  and  evaluation  by  comparing  their  theoretical  predictions  with  the  actual 
experimental  response  for  simple  as  well  as  complex  input  histories.  Some  of 
these  recent  comparative  investigations  are  those  of  Walker  /47/  (Walker  /38/, 

Krieg  et  al  /41/ ,  Miller  /44/),  Cemocky  /48 /  (Bodner  and  Partom/42/,  Cemocky 
and  Krempl  /37,  Lin  and  Wu  /35)  and  /49/  (Bodner  et  al  /42,43 /,  Cernocky  and 
Krempl  /37/,  Cernocky  /48/) ,  Beek,  Allen  and  Milly  /50/  (Allen  and  Haisler  /28 /, 
Valanis  :/.34/ ,  Cemocky  and  Krempl  /37/),  Walker  /38/,  Hart  /40/,  Krieg  et  al  /41/, 
Bodner  et  al  /42,43/,  Miller  /44,  Robinson  / 45/  and  Cescotto  and  Leckie  / 46/), 
and  Eftis  and  Jones  / 19/  (most  of  the  theories  in  / 26/-/46/).  Criteria  used 
for  evaluating  each  theory  include,  (a)  the  facility  of  its  calibration,  i.e. 
determination  of  material  constants  from  experimental  data;  (b)  its  predictive 
capabilities  of:  monotonic,  cyclic  hysteresis,  creep,  stress  relaxation  and 
multiaxial  responses;  (c)  its  sensitivity  to  numerical  integration;  and  (d)  the 
sensitivity  of  the  response  to  small  changes  in  material  constants. 

The  conclusions  of  the  comparative  study  of  various  constitutive  theories 
made  by  Beek,  Allen  and  Milly  / 50/  are  particularly  interesting.  These  are 
given  in  Tables  X,  XI  and  XII  taken  from  their  report.  The  uniaxial  formulation 
of  each  model  is  listed  in  Table  X,  using  a  common  notation.  Beek  et  al  / 50/ 
showed  that,  within  a  thermodynamic  framework,  all  models  are  consistent  with 
the  internal  state  variable  formulation.  Thus, any  model  can  be  written  in  the 
following  uniaxial  form: 

o  -  e  (e  -  -  gT)  (9) 

T 

where  °,  E,  g,  and  e  are  the  stress.  Young's  modulus, strain,  inelastic 
strain  (considered  as  an  internal  state  variable)  and  thermal  strain,  respectively. 
In  each  theory,  a  growth  law  for  the  internal  state  variable  ctj  is  further 
postulated,  and  of.  Is  obtained  by  time  integration.  Thus: 


-  26  - 


doj 

"dT*4!  (£>  T>  °2’  “3— •  “»>  <l0) 

and  *  .  »  » 

a l  m  J  a  (t  )  dt  ,  (11) 

*■ 

where  t  is  the  time,  T  :  the  temperature,  a ^  :  the  back  stress  which  is 
related  to  dislocation  arrangement  and  thus  produces  kinematic  hardening 
and  Bauschinger  effect,  and  :  the  drag  stress  which  represents  the 
dislocation  density  and  produces  isotropic  hardening. 

Differences  between  various  models  are  due  to  the  number  of  internal 
state  variables  and  their  growth  laws.  Also,  the  experimental  data  needed 
for  characterization  vary  greatly  between  models.  From  the  assessment  made 
by  Beek  et  al  /50/,  the  model  of  Cemocky  and  Krempl  /37/  seems  to  be  complex 
enough  to  model  viscoplastic  behavior  in  metals,  for  low  temperature  and  mono¬ 
tonic  load  conditions.  The  theory  proposed  by  Bodner  et  al  /42 ,43/  seems  to 
be  extremely  useful,  simple  to  construct  and  to  give  good  results  for  high 
strain  rate  applications  as  verified  recently  by  Bodner  and  coworkers.  This 
theory  has  also  been  applied  recently  by  the  present  author  with  satisfactory 
results  on  quasi-static  to  dynamic  strain-rate  change  in  tension  on  AISI  316 
stainless  steel  at  ambient  temperature  / 51/.  On  the  other  hand,  Stouffer  and 
Bodner  /52/  and  more  recently  Bodner  /53/  attempted  to  introduce  nisotropic 
hardening  into  the  original  formulation  of  the  theory.  Experimental  verifica¬ 
tions  of  these  modifications  have  not  been  reported  yet  in  the  literature. 

Based  on  the  previous  review  and  discussions,  it  may  be  concluded  that 
possible  models  to  recommend  for  fitting  the  present  experimental  results 
would  be  those  of  Cernocky  and  Krempl,  and  Bodner  and  associates.  Neither  of 
the  two  should  be  expected  to  give  complete  satisfactory  results,  however, 
mainly  because  of  the  large  plastic  strains  achieved  in  the  present  tests,  and 
accompanying  non-linearities  which  would  require  a  rigorous  theoretical  approach. 


SECTION  4 


CONCLUSIONS 


A  major  part  of  this  work  was  to  establish  new  data  bases  for  OFHC 
copper  and  plain  carbon  steel,  related  to  their  strain-history  dependence, 
particularly  in  terms  of  the  hysteresis  loop  upon  load  reversals  and  the 
associated  Bauschinger  effect,  at  quasi-static  as  well  as  impact  strain 
rates.  Such  information  should  be  accounted  for  and  used  to  remedy  the 
present  practice  in  which  many  fabricated  components  used  in  mechanical, 
civil  or  aeronautical  applications,  experience  loads  during  their  service 
life  which  pay  no  heed  to  previous  deformation  history. 

The  phenomenological  behavior  of  the  two  materials  was  obtained 
through  a  series  of  sequential  reverse  torsional  shear  tests  at  impact 
strain  rates  using  a  modified  split  Hopkinson  bar  and  a  quasi-static  strain 
race  using  a  slow  rotational  drive  unit.  The  observed  response  was  next 
evaluated  in  terms  of  various  Bauschinger  effect  parameters,  and  also  in 
the  light  of  predictions  from  classical  plasticity  theories  as  first  and 
crude  approximation.  Finally,  the  excellent  reviews  available  in  the 
literature  on  modem  constitutive  equations  recently  proposed  for  metals 
were  examined  and  those  most  probably  suited  to  model  the  present  results 
were  identified. 

The  main  results  of  the  present  investigation  ax£  hereafter  emphasized 
again  as  main  conclusions: 

(a)  Cyclic  loading  of  materials  at  intermediate  or  large  constant  strain 
amplitudes  and  at  impact  strain  rates  can  be  achieved  with  the  minimum 
amount  of  difficulty  in  torsional  shear  using  a  properly-modified 
classical  stored-torque  split  Hopkinson  bar  apparatus.  The  value  of 
the  strain  amplitude  can  be  controlled  by  the  position  of  the  friction 
clamp  along  the  system.  The  strain  rate  achieved,  and  hence  also  the 
strain,  is  essentially  controlled  by  the  amount  of  stored  torque  and 
also  the  specimen  gage  length. 

(b)  On  ductile  materials  like  aluminum  and  copper,  sequential  reverse  impact 
torsion  can  be  continued  for  a  large  number  of  cycles  until  instabilities 
limit  further  deformation.  With  materials  of  limited  ductility,  like 
steel,  this  ends  in  fracture  after  a  few  cycles. 
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(c)  Of  the  three  materials  tested,  the  aluminum  alloy  has  the  least  strain- 
rate  sensitivity,  under  monotonic  loading,  whereas  mild  steel  is  the 
most  sensitive.  This  is  in  line  with  our  current  knowledge  of  the 
strain-rate  sensitivity  of  metals. 

(d)  The  stress-strain  behavior  of  all  three  materials  is  affected  by  strain 
history.  The  effect  is  more  appreciable,  however,  for  highly  strain- 
rate  sensitive  materials. 

(e)  For  dopper,  irrespective  of  previous  loading  history,  the  flow  curve  in 
sequential  reverse  loading  assumes  a  lower  level  than  that  of  continuous 
forward  loading  starting  from  a  total  cumulative  strain  of  0.95.  In  mild 
steel,  however,  plastic  flow  following  reverse  yielding  proceeds  always 
at  a  level  below  that  expected  from  continuous  forward  loading.  For  both 
materials,  the  flow  curve  flattens  out  with  increasing  sequential  loading, 
ultimately  reaching  saturation, 

(f)  Reverse  loading  characteristics,  particularly  the  Bauschinger  phenomenon, 
can  be  described  at  all  strain  rates  by  means  of  quantitative  parameters, 
such  t  ,  Tg  and  t^/At.  In  general,  values  of  each  parameter  vary  with 
strain  rate.  Consequently,  it  is  inaccurate  to  describe  the  reverse 
loading  behavior  of  strain-rate  sensitive  metallic  materials  by  using  data 
collected  at  the  quasi-static  strain  rate,  as  it  is  often  followed  in  the 
analysis  of  structural  engineering  problems.  Instead,  it  will  be  much  more 
accurate  and  economical  to  use  values  of  the  yield  stress  ratio,  YSR  -  T  /x 
which,  for  a  given  prestrain,  do  not  vary  much  with  strain  rate  (cf.  Tables 
V  and  VIII). 


(g)  Comparing  the  present  experimental  results  with  predictions  of  the  ctfpssical 
anisotropic  hardening  theories  indicate  that  reasonable  agreement  does  not 
seem  to  exist  when  comparison  is  made  with  various  modified  versions  of  the 
original  Prager’s  kinematic  hardening  rule.  This  is  not  unexpected,  given 
the  fact  that  the  present  data  were  collected  at  fairly  large  strains,  high 
strain  rates  and  generally  non-isothermal  conditions. 

(h)  Reviewing  the  literature  on  advanced  constitutive  models  recently  proposed 

for  metals,  and  examined  in  more  details  by  other  investigators,  it  appears 

1* 


that  the  theories  of  Bodner  and  associates  and  of  Cernocky  and  Krempl  ofrer 
the  greatest  promises  for  fitting  the  present  experimental  data.  However, 
because  of  the  finite  strains  involved  in  the  present  data,  such  statement 


should  be  checked  by  actual  comparisons  In  order  to  provide  directions 
for  more  rigorous  theoretical  treatments,  If  needed.  This  should  be  a 
worthy  research  area  to  explore  in  future  studies  in  dynamic  plasticity. 
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TABLE  I  -  Composition  of  Materials  Tested 
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TABLE  V  -  Liot  of  Some  Parameters  Representative  of  Bauschinger  Effect  and  Hardening  for  Reverse  Loading  on  Copper 
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TABLE  VIII  -  List  of  Some  Parameters  Representative  of  Bauschinger  Effect  and  Hardening  for  Reverse  Loading 
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Represents  the  maximum  stress  level  reached  on  the  immediate  previous  loading 
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TABLE  XII  Required  Material  Parameter  Characterization  in  Various  Viscoplastic 
Constitutive  Equations  (Beek,  Allen  and  Milly  /50/) 


Cernocky  and  Krempl 

Constant  Strain  Rate  Tensile  Tests  with 
Intermittent  Hold  Times 

Relaxations  Tests 

Krieg,  Swearengen, 

Stress  Drop  Tests 

and  Rohde 

Constant  Strain  Rate  Tensile  Tests 

Bodner  et  al. 

Constant  Strain  Rate  Tensile  Tests 

Creep  Tests 

Walker 

Constant  Strain  Rate  Cyclic  Tests 

Constant  Strain  Rate  Tensile  Tests 

Miller 

Creep  Tests 

Constant  Strain  Rate  Cyclic  Tests 

Constant  Strain  Rate  Tensile  Test 

Cescotto  and  Leckie 

Constant  Strain  Rate  Cyclic  Tests 

Constant  Strain  Rate  Tensile  Tests 

Stress  Drop  Tests 

Hart 

Relaxation  Tests 

Robinson 

Stress  Drop  Tests 

Valanis 

Constant  Strain  Rate  Tensile  Test* 

Allen  &  Haisler 

Constant  Strain  Rate  Tensile  Tests  At 
Several  Constant  Temperatures 

Creep  Tests  At  Several  Constant 
Temperatures 

Single  Cycle  Strain  Controlled  Test 

At  Each  Constant  Temperature 

Schematic  drawing  of  torsional  SHB  with  insets  showing  various  details 
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Specimen 
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Camera 
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Fig.  2  Closeup  of  CRO  and  related  instrumentation 


Dimensions  in  mm 


Fig.  7  Configuration  and  nominal  dimensions  of  test 
specimens 
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Fig.  11  Typical  records  of  quasi-static  reverse  loading 
copper  specimens 


Fig.  12  Typical  records  of  quasi-static  reverse  loading 
steel  specimens 
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Fig.  14  Typical  oscilloscope  records  of  impact  reverse  loading  on  steel 
specimens 


under  sequential  reverse  loading  at 


Shear  Strain  ,  y 

Fig.  15(c)  Stress-strain  curves  for  copper  under  sequential  reverse  loading  at  i  .  1100 
constant  strain  amplitude  y  ~  0.97 


constant  strain  amplitude 


Shear  Strain  ,  y 

Fig.  16(b)  Stress-strain  curves  for  copper  under  sequential  reverse  loading  at  y  -  550  s"1,  for  a  constant 
strain  amplitude  y  ~  0.32 


Shear  Stress  ,  'T(MPa) 


Shear  Sirain ,  y 


Fig.  17 (b)  Stress-strain  curves  for  copper  under  forward 
and  reverse  loading  at  y  *  10~3  s”l,  for  a 
constant  strain  amplitude  y  -  0.76 


Shear  Strain  ,  y 

Fig.  17(c)  Stress-strain  curves  for  copper  r-ider  forward 
and  reverse  loading  at  Y  -  in~3  s~^-.  for  a _ 
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Fig.  19  Variation  of  reverse  yield  stress,  t  ,  and  Bauschinger 
stress,  Tg,  for  copper,  with  strain  Amplitude,  yq,  for 
various  strain  rates 
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Fig.  20  Variation  of  Bauschinger  effect  indicator,  tb/At,  for  copper 
with  strain  amplitude,  y  ,  for  various  strain  rates 
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PART  II 


Ratchetting  of  Pressurized  Tubes 
Subjected  to  Radial  Cyclic  Thermal  Gradient 


SECTION  1 


INTRODUCTION 


The  structural  ratchet  ting  phenomenon  ( cumulative  cyclic  plastic 
strain  growth),  being  of  prime  importance  in  modern  systems  where  a 
steady  and  a  cyclic  load  components  are  acting  simultaneously,  has  been 
investigated  within  the  frame  of  the  present  work.  This  phenomenon  and 
its  associated  ratchet  strains  are  greatly  affected  by  the  formulation 
of  the  structure,  the  flow  and  hardening  characteristics  of  its  material. 

If the  cyclic  load  component  is  due  to  a  cyclic  thermal  gradient,  ratchet 
strains  are  also  affected  by  the  form  of  that  gradient  as  well  as  its 
amplitude. 

In  the  particular  case  of  a  steel  tube  subjected  to  steady  internal 
pressure  and  cyclic  radial  thermal  gradient,  which  simulates  the  conditions 
of  fuel  cans  of  liquid  Metal  Fast  Breeder  Reactors,  ratchetting  -  among 
other  deformational  patterns  -  may  occur  under  certain  loading  combinations/1/. 

All  previous  analytical  work  on  this  problem  was  confined  to  Bree's 
uniaxial  formulation  of  a  tube  made  of  an  elastic  -  perfectly  plastic  or 
elastic  -  linearly  plastic  material  with  kinematic  hardening  assumptions. 

In  the  present  phase,  an  effort  was  made  to  provide  a  more  faithful  biaxial 
formulation  of  the  tube  problem  and  to  develop  a  new  numerical  procedure 
capable  of  handling  any  tube  dimensions,  different  temperature  distributions, 
and  various  hardening  models.  This  new  formulation  has  been  presented  in  a 
previous  report/2/,  and  will  only  be  summarized  hereafter.  More  emphasis 
will  be  made  however,  on  the  numerical  procedure  developed  for  Its  implemen¬ 
tation  and  on  the  results  obtained  as  contrasted  with  those  of  the  uniaxial 
model  and  those  obtained  from  experiments. 
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SECTION  2 

BIAXIAL  FORMULATION  OF  THE  TUBE  PROBLEM 

The  elasto-plastic  investigations  to  determine  stresses,  strains 
and  modes  of  deformations  for  tubes  subjected  to  both  steady  internal 
pressure  and  cyclic  radial  temperature  gradient  as  previously  done  by 
Miller/3/.,  Bree  /4/,  Mulcahy  / 5/  and  others  were  confined  tc  thin-walled 
tubes  for  which  an  equivalent  uniaxial  model  was  suggested.  The 
hoop  stress  was  maintained,  the  thermal  stress  was  magnified  by  a  factor 

equal  to  ^y~y  and  the  axial  stress  was  neglected.  The  arguments  made 
by  Bree  showed  that  such  an  equivalent  uniaxial  model  -  for  the  originally 
two  dimensional  problem  -  would  lead  to  more  conservative  results,  since 
the  effect  of  the  axial  stress  -  which  was  neglected  -  is  opposing  that  of 
the  hoop  stress  and  therefore  the  results  obtained  were  expected  to  be 
greater  than  what  would  be  obtained  by  an  exact  solution. 

The  actual  elasto-plastic  solution,  being  in  two  dimensions,  is 

analytically  impossible.  Therefore,  from  experienc  -  earned  from  numerical 

investigations  of  the  uniaxial  formulation  of  the  tube  problem  / 1/  it  is 

realized  that  a  similar  procedure  could  be  followed  to  determine  the  elasto- 

plastic  response  of  thick  tubes  under  generalized  plane  strain  conditions. 

In  this  case,  the  radial  stress  o^  is  no  longer  neglected,  the  axial  stress 

is  not  set  to  zero,  the  thermal  stress  is  not  magnified  and  the  axial 

strain  e  is  assumed  to  be  constant, 
z 

A  complete  numerical  procedure  is  developed  and  presented  in  this  work 
for  the  generalized  plane  strain  tube  problem.  Different  hardening  rules 
for  the  tube  material  can  be  employed  in  this  numerical  procedure.  Both  the 
kinematic  and  isotropic  hardening  models  showed  satisfactory  performance  as 
related  to  the  tube  problem. 


2.1  .  iera]ized  plane  strain  formulation 

Consider  a  thick-walled  cylinder  with  inner  radius  a,  outer  radius  b 

subjected  to  steady  internal  pressure  p  and  a  cyclic  temperature  gradient 

T(r).  The  tube  has  closed  ends  and  is  further  subjected  to  an  additional 

uniform  axial  stress  oq.  Generalized  plane  strain  assumptions  imply  that 

the  axial  strain  e  is  independent  of  radius  r,  i.e.  e  =  constant.  Tota 
z  z 

strain  is  composed  from  elastic,  plastic  and  thermal  components.  Creep 
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As  shown  in  detail  elsewhere  HI ,  equilibrium,  compatibility  and 
stress-strain  relations  can  be  used  to  derive  the  following  expressions 
for  the  stress  components: 


a2  b2-  r2 

°r  =  -P*F-*^T^r  + 


fl  +  I  +  I  -  fa2  r2_a2  T  1 
U1  l2  l2  7T'W^ar  U' 


a2  b2+r2  E 

°6  "  p  r2  *b2-a2  (l+v)(l-2v) 


-'2-'3-7r 


b2  r2+a2 


)2-a2  V 


•  K  -  v(e?+eo)  +  U(1+V)]T  , 


°z  =  2vPPTF  + 


Wbr  lo>- T^?Lc0-v(^+GS)] 


and  strain  components: 


(l+v)(l-2v)  i  r-.  b2  1  1 1  i  i  !?+I3  ^ 

"Efbi-aaj  ”pa  +  T^T  “  F^iT 


TTl^Wl  '  VGz  +  T^J  ^1_2v)  er  +  a  (1+v)  T1 

(l+v)(l-2v)  2  r»  ,  b2  1  1  .  I2+I3  b2lo 

f.(b2-a2)  p  a  P  +  F  +  l)  ‘  TTv  '  b'~ 


t1  +hj^  -  VV 


(2) 
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where 


r  eP  .  £P 

j  -  1-2v  (  cr  C6 


‘l  '  IfT^I 


/ 


dr. 


‘2  *Xf$F2  1  r  {€r*  e8>  dr- 


u 

I.  .  -c.(l-2yj(Hv)  ;rTrdr 
3  (l-v)r2  , 


(3) 


and 


’  <‘l  +  ‘2  *  '3!^ 


Note  that  the  value  of  e 
This  constant  value  of  e 
equilibrium,  as: 


z 

z 


is  an  unknown  constant  yet  to  be  determined, 
is  determined  from  the  condition  of  axial 


cz  = 


T 


(l-2v) 


pa1  .  _ 2v _ 

(b*-a‘)£  ( \+v){ l-2v) (bl-al) 


b 

[b2IQ-2  /  rl-jdr] 

d 


(4) 


In  the  above  formulation,  the  uniform  axial  stress  Oq  can  be  set 
to  any  value  and  hence  this  feature  allows  investigating  plane  strain 
conditions  (e^  =0)  since  the  corresponding  value  of  a ^  can  be  determined 
from  Eq.  (4) 

2. 2  Plasticity  relations 

Since  the  tube  is  subjected  to  cyclic  thermal  gradient,  the  material 
model  to  be  used  in  the  analysis  should  exhibit  the  characteristics  of 
metals  when  subjected  to  cyclic  stresses.  The  most  important  of  these 
characteristics  for  initially  annealed  metals  and  alloys  are  cyclic 
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hardening,  cyclic  creep  and  cyclic  relaxation.  As  a  first  step  toward 
the  modelling  of  these  features  of  material  behaviour,  the  tube  is 
analysed  here  using  classical  plasticity  theories  with  both  the  kinematic 
and  isotropic  hardening  rules. 

2.2.1  Kinematic  hardening 

The  plasticity  relations  for  kinematic  hardening  can  be  formulated 
as  follows.  The  yield  condition  for  axisymmetnc  states  of  stress  is 
given  by: 

F  =  ($r-a.)’  ♦  (S^)2  ♦  (Sz-a2)2  .  §  =  0  (5a) 


were  S^,  and  are  the  stress  deviators,  given  by 

sr  =  or  -  <<v«e«2>  /  a  ,  (ft) 

se  *  °e  -  <VW  1  3  • 

S2  =  o2  -  (or«e«z)  /  3  , 

and  a  ,  a  ,  a  are  the  coordinates  of  the  center  of  the  yield  surface 

IT  8  z  t 

in  stress  space.  The  initial  yield  stress  is  denoted  by  q.  The 
flow  rule  -  which  relates  plastic  strain  increments  to  stresses  -  is 
based  upon  Drucker  postulate  which  states  that  the  plastic  strain  rate 
vector  is  normal  to  the  yield  surface,  so  that  : 


dep  =  dX 
r 


dep  =  dX 


3F  _  2 
3Tr  -  3  dX 

3F  .  2 
3oe ' 3 


.  p  j,  3F  2 
de;  =  dX  .  T  dX 
z  3a  3 


(Sr  -  ar)  . 


(se  "  °e)  » 


(s2  -  ,z)  , 


where  dX  is  a  positi”.-.  nlar  multiplier  yet  to  be  determined.  Introducing 
an  effective  plastic  &  rain  increment  deP  which  is  defined  as  : 


de?  =  (-1  (de°  .  dep  -  d£:p  .  dep  -  de;  .  de;)] 


(7) 
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‘,/j 
V.'J 
v;j 
V  'j 


If 


which  reduces  de^  to  the  observed  plastic  strain  in  the  uniaxial  test  and 

3  2  dcp 

substitution  of  Eq.  (6)  in  Eq.  (7)  yields  the  value  of  dX  as  (•?)  — r-  • 

1  o 

Hence,  the  flow  relations  can  be  rewritten  as: 


der  -  7 


dr 


a 

& 


(Wl 


P  .  1  IE!  fs  .o  1 

a  ~  i  •  _  1^6  o' 


de 


dcz  =  I 


dcp 


(8) 


The  hardening  rule  whi  \  describes  the  incremental  translation  of  the 
center  of  the  yield  surface  is  based  upon  Prager  kinematic  hardening 
model  which  assumes  that  such  translations  are  in  the  direction  of 
the  plastic  strain  increments-  Thus: 


>2 


L- 

r. 


da  =  du  .  de;,  ; 

r  r  * 


da0  =  du  -  ; 


da,  =  du  .  dci 


(9) 


where  du  is  a  positive  scalar  multiplier.  Introducing  an  effective 
incremental  translation  da  of  the  center  of  the  yield  surface  as 


da  =  [|  ,  (d«r  .  da^  +  daQ  .  daQ  +  daz  .  daz)]5 


(10) 


reduces  da  to  the  observed  incremental  translation  of  the  center  of 

the  yield  surface  in  the  uniaxial  test.  Substitution  of  Eq.  (9)  into 

do 

Eq.  (10)  yields  dy  =  - which  in  turn  gives  da^,  daQ  and  daz,  as 

d  sP  r 


da  .  1  .  -22. 
r  3  dlp 

„„  Z  da 

%-7-^p 

da  =  §  .  4 
2  3  dcp 


dep  =  |  C  .  dep  . 


dep  =  |  C  .  dtp  . 
dcp  =  |  C  .  ddp  . 


(ID 


.V.* 
4V 
.v v 

k 


.V 

*  V 
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where  C  = 


is  the  coefficient  of  strain  hardening  as  determined 


from  the  basic  uniaxial  test. 

The  yield  condition,  Eq.  (5),  the  flow  rule,  Eq.  (8),  and  the  hardening 
rule,  Eq.  (11),  completely  define  the  plasticity  theory  for  kinematic  harden¬ 
ing.  In  the  present  work,  the  coefficient  of  strain  hardening  C  is  taken  as 
a  constant  which  implies  that  linear  hardening  is  assumed. 

2.2.2  Isotropic  hardening 

If  the  isotropic  hardening  is  considered,  the  yield  condition,  Eq.  (5), 
becomes 

F  =  Sr  +  s0  +  Sz  -  I  ay  =  0  (12) 

where,  in  this  case,  o^  is  the  current  value  of  the  yield  stress  which  is 

equal  to  the  last  value  of  the  equivalent  stress  o  at  which  plasticity  has 

occurred.  The  initial  value  for  o  is  the  initial  yield  stress  o 

y  J  y,o 

The  flow  rule  of  Eq.  (8)  becomes 


3  de*3 


de0  "  1  ““  se 
a 


d£p  =  1  „ 

z  2  a 


The  hardening  rule  in  the  case  of  isotropic  hardening  is  simple,  since  no 
translation  of  the  center  of  the  yield  surface  occurs,  i.e.  =  oQ  =  =  0, 

and  the  current  value  of  the  yield  stress  is  set  equal  to  the  value  of  the 
equivalent  stress  a  at  which  plasticity  has  occurred  for  the  last  time. 

The  yield  condition,  the  flow  rule  and  the  hardening  rule  are  employed 
in  the  numerical  procedure  for  either  of  the  kinematic  and  isotropic  material 
models  and  the  solution  for  stresses  and  strains  is  then  found  at  any  point 
in  the  tube  during  any  given  cycle. 
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SECTION  3 

NUMERICAL  PROCEDURE 


Cyclic  thermo-elasto-plastic  problems  may  be  solved  by  means  of  the 

method  of  successive  elastic  solutions  in  which  the  total  thermal  gradient 

is  applied  incrementally.  For  each  increment,  an  iterative  procedure  is 

employed  in  order  to  determine  the  plastic  strain  rate  vector  (de*3,  dc*3,  dep) 

-  .  r  0  z 

for  the  biaxial  tube  model  which  satisfies  the  equilibrium,  compatibility 
and  stress-strain  relations.  Then,  this  plastic  strain  rate  vector  is  used 
to  determine  the  corresponding  stresses. 

The  procedure  used  for  the  uniaxial  tube  model  is  basically  the  same 
as  that  used  for  the  biaxial  tube  model.  Therefore,  the  numerical  procedure 
for  the  uniaxial  case  (being  a  special  case  of  the  biaxial  model)  is  not 
given  here  and,  instead,  that  for  the  biaxial  problem  is  presented  in  full. 

As  mentioned  earlier,  the  linear  temperature  distribution  is  only  a 
special  case  of  the  bilinear  distribution.  Assigning  zero  to  the  value  of  'h' 
in  Fig.  1  will  convert  the  bilinear  temperature  distribution  into  a  linear 
distribution.  Hence  the  proper  value  of  Th’  within  the  following  numerical 
procedures  will  provide  the  different  versions  for  the  different  temperature 
distributions. 


The  procedure  may  be  divided  into  two  main  segments;  (a)  determination 
of  the  temperature  level  at  which  first  yield  occurs,  and  (b)  continuation 
of  thermal  load  incrementation  to  determine  the  elasto-plastic  solution  of 
the  problem. 


3. 1  Determination  of  the  Temperature  Level  at  Which  First  Yield  Occurs 

For  the  general  case  in  which  a  bilinear  temperature  distribution  is 
considered  across  the  tube  wall,  the  temperature  may  be  given  by 


T(x)  =  0 
T(x)  = 


o  <  x  h 

h  <  x  <  1 


(14) 


where  x  = 


r)  ft. 


-  T(x-h)  /  (1-h) 
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The  temperature  changes  proportionally  at  all  radii  such  that,  for  an 
intermediate  temperature  as  shown  by  a  dotted  line  in  Fig.  1,  AT  is  replaced 
by  AT'.  For  the  full  temperature  distribution  the  temperature  difference  is 
AT. 


In  order  to  determine  the  temperature  level  at  which  first  yield  occurs, 
the  total  temperature  gradient  AT  should  be  divided  into  an  appropriate  number 
of  small ‘increments.  Care  should  be  given  in  choosing  the  number  of  load 
increments  to  assure  convergence  of  the  iterations  and  to  save  computer  time. 
According  to  Barsoum  /6/,  and  also  from  numerical  experiences,  it  was  found 
that  relatively  large  load  increments  would  cause  divergence.  The  reason 
behind  this  divergence  is  that  the  corrected  plastic  strains  are  calculated 
through  stresses  which  are  in  turn  erroneous  during  the  iteration  process. 

The  thermal  load  is  now  applied  incrementally  and  the  thermo-elasto- 
plastic  problem  is  solved  with  the  initial  values  of  plastic  strains 
"r  =  =  £z  =  0  (°r  ^eP  “  0  for  uniaxial  model).  After  each  temperature 

increment,  a  check  is  made  on  the  yield  function  F(r)  at  different  points 
covering  the  thickness  from  the  inner  to  the  outer  radius.  If  F(r)  <  0,  no 
yield  has  started  yet.  If  F(r)  >_  0,  yield  has  occurred;  consequently,  the 
iterative  procedure  for  evaluating  the  plastic  strains  is  to  be  started,  and 
the  thermal  load  incrementation  is  Cv..tinued  until  full  temperature  difference, 
AT,  is  cover. 

3.2  Continuation  of  Thermal  Load  Incrementation  to  Determine  the  Elastic 
Plastic  Solution 

For  the  evaluation  of  correct  stresses  and  strains,  the  procedure  used 
consisted  of  the  following  steps. 

(1)  For  all  radii  at  which  plastic  deformation  is  likely  to  occur,  initial 
values  for  the  increments  of  plastic  strains  are  assumed  such  that 
Ae^+Atg  +  Ae^  =  0.  An  appropriate  starting  assumption  is  to  take  these 
increments  very  small,  as  for  instance 

AeP  =  1.5  x  lCf6,  Ae?  =  -  1  x  10'6,  AeP  =-0.5  x  10“6  (15) 

r  u  2 


f  %, 

u 


v^* 


«/• 

•V 

v: 


.  * « 
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It  should  be  noted  that  Ac**,  AeP  and  AeP  are  taken  as  zeros  at  all 

r*  0  z 

other  radii  at  which  no  yield  occurs. 

(2)  At  every  radius,  the  effective  increment  of  plastic  strain,  As*5,  is 
calculated  as 

1 

AeP  =  [|  (Aep2  +  A£p2  +  AEp2)]*  (7)  bis 

(3)  At  every  radius,  the  change  in  the  coordinates  of  the  center  of  the 
yield  surface  is  calculated  as 

=  f  C  *  Aer  ;  =  |  C  .  Aep  ;  Actz  *  |  C  .  Aep ,  (11)  bis 

BE 

where  C  is  the  work-hardening  parameter  taken  as  -r—r 

1  — P 

(A)  At  all  radii,  the  new  values  of  plastic  strains  are  calculated  as 


■P 

r,n 

=  er,n-l 

+ 

AeP  , 
r  * 

.P 

9,n 

=  e6,n-l 

+ 

■ 

(16) 

P 

z,n 

=  ez,n-l 

+ 

4eS  ’ 

where  n-1  and  n  denote  a  previous  and  current  increment  of  thermal  .  \ 

loading. 


(5)  At  all  radii,  the  new  values  of  the  coordinates  of  the  center  of 
the  yield  surface  are  calculated  as 


a  = 

°r,n-l 

+ 

Act  , 

r»n 

afi  n  = 

o,n 

ae,n-l 

+ 

A0e  . 

(17) 

a  = 

az,n-l 

+ 

Act 

z.n 

z 

(6)  Using  »  ef  _  and 

r,n  0,n  z 

obtained 

»n 

in  step  4, 

the  elasto-plastic 

thermal  problem  is  solved* using  Eqs.  (1),  (2)  and  (4)  to  calculate 

the  current  stresses  and  strains:  a  ,  o„  ,  e  -  e  and  £  .At 

r,n  0,n  z,n  r,n’  9,n  z,n 

this  level  of  calculations,  it  should  be  noted  that  the  integrals  and 
1^  will  have  values  other  than  zero. 


(7)  At  all  radii,  the  deviatoric  stress  tensors  are  calculated  as  : 

_  a.  +  a 

~  _  2  r.  Q,n  z,ni 

s„  .  *  1°  n - 1 - 5 — H  ’ 


!  Gv.n 

!  [°9,n 


!  [«z.n 


0  +0 

JVL_LJ1], 


a.  +  o 
6,n  r,ni 


(5b)  bis 


Using  the  values  of  ex  ,  a*  and  a  obtained  in  step  *5* 
6  r,n  0,n  z,n  r 

together  with  S  ,  S  and  S  obtained  in  step  ’7',  the 
°  r,n*  0,n  z,n  f  * 

A  P*  A  P* 

corrected  values  of  plastic  strain  increments  Aer  ,  Ae!:  and 
P* 

Aer  are  as  follows  161  : 


,_P*  =  E  •  a£P*  -  E  •  a£P*  =  ^  e  (18) 

r  e  *  erl’  Ae0  *  01  *  Aez  eo  *  zl (18) 

e  e  e 


where  is  an  equivalent  modified  total  strain  calculated  as 

Ee  =  f|  <cri  +  e61  *  ezl>lJ  (19) 

and  Ae11  in  this  case  is  calculated  from 


-p  2  1+v  _ 

Ac  *  ee  '  3  '  T  °y,o 


where  a  is  the  initial  yield  stress,  while  e  ,,  and  e  , 

y,o  J  *  rl’  01  zl 

are  the  deviatoric  modified  strains  to  be  calculated  from 


rl 

1+V  r 
=  1 

•Sr,n 

ar,n^ 

+ 
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(9)  The  corrected  value  of  effective  plastic  strain  increment  is 
calculated  at  all  radii  using  the  new  values  of  plastic  strain 
increments  obtained  in  step  '8* 

tep‘  =  [f  t(4ep*)2  +  (4epV  +  (icp*)2)]!  (22) 

(10)  The  convergence  criterion  is  checked  at  all  radii  : 

|AeP*  -  AeP|  <  n  (23) 

where  n  is  a  constant  to  be  chosen  on  basis  of  accuracy  required 
and  availability  of  computational  time,  since  a  smaller  n  will 
increase  the  number  of  iterations  performed  to  reach  the  solution 
and  hence  computational  time  is  increased.  AeP  in  the  above  conver¬ 
gence  criterion  is  the  value  obtained  in  step  *2'.  If  the  conver¬ 
gence  criterion  is  satisfied  at  all  radii  at  which  yield  is  occuring 
proceed  to  step  'll',  and  if  the  convergence  criterion  is  not  satisfied 
go  to  step  *2*  and  repeat  the  calculation  (step  '2'  to  step  *10*) 

using  the  newly  calculated  values  of  plastic  strain  increments,  viz, 

P*  P*  .  P* 

Ae  ,  Ae.  and  Ae  obtained  in  step  8  until  convergence  is 

T  u  Z 

attained. 

(11)  At  all  radii  which  were  not  yielding  during  the  thermal  load  increment 
considered,  the  yield  function  F(r)  is  calculated  from  Eq.  (5a)  to 
determine  whether  yield  has  occurred  during  the  current  thermal  load 
increment.  If  yield  occurs,  these  radii  should  be  allowed  to  deform 
plastically  during  the  application  of  subsequent  thermal  load  increments. 

(12)  The  previous  steps  are  repeated  until  the  full  temperature  distribution 
is  applied  to  the  tube. 

In  the  above  numerical  procedure,  the  relations  for  a  kinematic  hardening 
material  were  employed.  If  any  other  hardening  model  is  used,  exactly  the  same 
procedure  may  be  applied  with  proper  employment  of  the.  hardening  rule,  flow 
rule  and  yield  function. 


As  regards  the  isotropic  hardening  model  for  the  material,  the  above 
procedure  is  followed  with  minor  modifications.  Steps  '3'  and  '5'  are 
deleted.  In  step  ’8*  the  corrected  values  of  plastic  strain  increments 
are  calculated  as  follows  /!/: 


CP*  _  Aep  .  p*  _  Aip  •  eP*  -  _  (iQ)  bis 

er  ~T  ’  rl  *  e6  *  "7"  *  £ei  ’  Z  *  T“  ‘  ezl  ’  (  )  biS 

.  e  ^  e 


where  ee  is  the  equivalent  modified  total  strain  calculated  as  : 


:e  =  If  <erl!  *  eei*  *  hi  >1* 


(19)  bis 


and  Ae^  is  calculated  as  : 


AiP  * 


.  [4cp  .  |  (^)(5  .  0y)] 


AeF  :  is  as  calculated  in  step  *2’ 

°  :  is  the  equivalent  stress  at  the  current  point, 

°y  :  is  the  current  yield  stress  at  the  current  point, 

and  ee^»  e@1  and  e  are  the  deviatoric  modified  strains  calculated  as 

en  *  (Sr,n}  +  Ae?  * 


e91  =  ^e.n5  +  Ae«  ’ 


ezl  "  T*  <Sz,n>  +  A£Pz  . 


In  step  ’ll',  the  yield  function  will  be  as  given  by  Eq.  (12;  and 
in  addition  to  checking  points  that  were  not  yielding  during  the  load 
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increment  considered,  the  yield  stress  of  the  material  should  be  corrected 
for  points  at  which  yielding  has  occurred  and  convergence  is  reached.  For 
those  points,  and  before  proceeding  to  the  next  thermal  load  increment,  the 
new  yield  stress  is  calculated  as 


For  other  points  where  yield  has  not  occurred  during  this  load  increment, 
the  yield  stress  should  remain  unchanged.  These  steps  are  repeated  until 
the  full  temperature  distribution  is  applied  to  the  tube. 

The  flow  chart  shoxm  in  Fig.  2,  indicates  the  main  portions  and 
features  for  the  biaxial  program  under  kinematic  hardening  assumption. 

For  isotropic  hardening,  the  main  features  are  basically  the  sane  except 
for  the  relations  according  to  which  the  new  increments  of  plastic  strains 
are  calculated. 


SECTION  4 

RESULTS  AND  DISCUSSIONS 

4.1  Choice  of  Numerical  Experiments 

The  two  computer  programs  developed  for  handling  the  tube  problem, 
based  on  the  uniaxial  and  biaxial  models  respectively,  were  carefully  tested 
and  tried  for  load  sets  causing  different  regimes  of  deformation.  The  two 
programs  showed  successful  converging  results  for  the  kinematic  and  isotropic 
versions  of  each  program. 

A  set  of  actual  material  constants  were  fed  into  these  programs.  These 
were  chosen  for  the  304  SS  used  by  Corum  and  Sartory  in  their  tube  ratchetting 
tests  18/ .  Thus,  comparison  between  numerical  and  experimental  results  was 
made  possible.  The  properties  of  304  SS  at  an  elevated  temperature  of  900°F 
(482°C)  are  represented  by  the  following  quantities  : 

,  g  CO 

Young’s  modulus,  E  =  23.3  x  10  psi  (16.38  x  10  kp/cm4') 

Poisson’s  ratio,  v  =  0.3 
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Yield  stress,  a *  11.15  x  10^  psi  (783.8  kp/cm2) 
Coefficient  of  strain  hardening,  8  *  0.0253 

Coefficient  of  linear  thermal  expansion, a  =  11.09  x  10  ^  °F  ^(19. 96x10  ^  °C  * 

Tube  dimensions  representing  Corum’s  test  specimen  were  also  used  in 
these  programs.  These  were  as  follows  : 

Inner  radius,  R^  =  3.845  in.  (97.663  mm) 

Outer  radius,  R  =  4.22  in.  (107.188  mm) 

o 

Mean  radius  to  thickness  ratio,  R  /t  *  10.75 

m 

In  order  to  complete  the  information  needed  to  run  either  one  of  the 
two  programs,  the  values  of  the  internal  pressure  p,  the  maximum  temperature 
difference  AT,  the  number  of  steps  by  which  this  difference  is  to  be  incremen¬ 
ted,  and  the  value  of  ’h’  which  defines  the  point  at  which  the  two  straight 
lines  of  the  bilinear  temperature  distribution  intersect  (h  =  0  for  linear 
distribution),  were  also  needed. 

The  strategy  followed  to  get  the  numerical  results  was  as  follows  : 

(a)  Initial  investigation  of  the  validity  of  the  adopted  numerical 
procedure.  The  cornerstone  for  this  validation  was  the  kinematic 
hardening  version  of  the  uniaxial  program.  A  complete  study  was 
made  on  this  version,  since  it  is  the  only  one  to  which  a  correspond¬ 
ing  analytical  solution  is  available  /9/.  Different  sets  of  loading 
parameters,  as  shown  in  Table  I,  which  covered  the  full  interaction 
diagram  for  the  specific  tube  considered.  Fig.  3,  were  used. 

(b)  Upon  verification  of  the  validity  of  the  kinematic  hardening  version 
of  the  uniaxial  program,  the  segment  of  the  program  concerning  the 
hardening  rule  was  modified  to  accommodate  the  isotropic  hardening 
rule.  Some  loading  samples,  as  shown  in  Table  II,  were  then  tried. 

(c)  Then,  being  sure  that  the  adopted  numerical  procedure  was  in  good 
order,  the  uniaxial  program  was  further  modified  to  yield  the  biaxial 
versions  for  both  kinematic  and  isotropic  hardening  models,  with  the 
numerical  technique  maintained  unchanged. 


(d)  The  biaxial  versions  were  then  used  to  get  results  for  the  pressure 
and  temperature  difference  used  by  Corum  in  his  test  (p  *  700  psi, 

AT  -  167°F). 

Several  test  runs  were  done  for  each  of  the  versions  used  until 
satisfactory  performance  was  achieved  and  enough  experience  was  earned. 

Only  computer  runs  of  interest  were  completed  and  printed.  Results  are 
presented  in  the  following. 

4.2  Results 

For  each  of  the  loading  sets  mentioned  above  the  corresponding  version 
of  the  program  was  used  to  calculate  stresses  and  strains  at  51  equally- 
spaced  points,  covering  the  full  thickness  of  the  tube  for  the  uniaxial 
versions,  and  26  points  for  the  biaxial  version.  The  resulting  stresses 
and  strains  were  printed  after  each  half  thermal  cycle  and  for  30  complete 
cycles . 

4.2.1  Kinematic  hardening  results 

(a)  For  a  pressure  p  of  100  psi  and  a  temperature  difference  AT  of  50°F 

at  the  inner  tube  surface  (x»l).  Fig.  4  shows  the  resulting  stress  distribu¬ 
tions  for  the  first  cycle  (solid  for  the  cooling  half  cycle  and  dotted  for 
the  heating  half),  and  the  growth  of  total  strain  e  with  cycling.  During  the 
cooling  half  of  the  first  cycle,  plastic  deformation  occurs  only  at  the 
inner  side  of  the  tube  (from  x  *  0.92  to  x  *  1),  and  no  plastic  deformation 
is  observed  during  the  second  half  of  the  first  cycle.  Further  cycling  shows 
no  growth  of  plastic  strains  and  the  stress  distributions  are  exactly  the  same 
as  in  the  first  cycle.  The  cyclic  steady  state,  as  could  be  seen  from 
Fig.  4(a),  (b),  is  then  an  elastic  shakedown.  This  agrees  also  with  analytical 
predictions  for  the  pressure  and  temperature  difference  employed  (as  could  be 
seen  from  Fig.  3).  The  value  of  the  plastic  front  (x  =  0.92),  as  well  as 
stresses  and  strains  at  all  points  are  also  found  to  be  the  same  as  calculated 
theoretically. 

(b)  For  p  »  100  psi  and  AT  »  200°F.  Fig.  5  gives  the  resulting  stresses  for 
the  first  cycle  and  the  growth  of  the  total  plastic  strain  with  further  cycling. 
During  the  cooling  half  of  the  first  cycle,  compressive  plastic  deformations 
take  place  at  the  outer  tube  side  from  x  =  0  to  x  =  0.36,  meanwhile  tensile 


plastic  deformations  occur  at  the  inner  side  from  x  *  0.58  to  x  =  1.  During 
the  heating  half  of  this  cycle  only  compressive  plastic  deformations  take 
place  at  the  inner  tube  side  from  x  *  0.72  to  x  =  1.  Further  cycling  gives 
exactly  the  same  stress  distributions  as  in  the  first  cycle.  Hence,  the 
range  from  x  =  0.72  to  x  *  1  is  undergoing  cyclic  plasticity  and  no  total 
strain  growth  is  observed.  The  cyclic  steady  state  is  thus  cyclic  plasticity 
F, ,  and  has  been  reached  after  the  first  cycle.  The  analytical  results  for 
the  same  operating  conditions  and  corresponding  regime  of  deformation  were 
found  to  be  exactly  the  same  as  the  numerical  predictions. 

(c)  For  p  =  100  psi  and  AT  »  300°F.  Results  are  shown  in  Fig.  6.  During 
the  first  half  cycle,  plastic  strains  are  developed  at  both  sides  of  the 
tube;  from  x  =  0  to  x  =  0.4  compressive  yielding  occurs,  and  from  x  =  0.56  to 
x  *  1  tensile  yielding  takes  place.  During  the  second  half  cycle  plasticity 
is  still  seen  at  both  sides;  tensile  from  x  =  0  to  x  *  0.38  and  compressive 
from  x  -  0.65  to  x  *  1.  Further  cycling  gives  exactly  the  same  stress  dis¬ 
tributions  seen  in  the  first  cycle.  Hence  the  regions  from  x  «  0  to  x  *  0.38 
and  from  x  =  n.66  to  x  =  1  are  undergoing  cyclic  plasticity  while  no  total 
strain  growth  is  observed.  The  cyclic  steady  state  in  this  case  is  cyclic 
plasticity  F^  and  is  reached  after  the  first  cycle.  These  predictions  agree 
also  with  analytical  predictions  for  the  same  operating  conditions. 

(d)  For  p  -  100  psi  and  AT  *  100°F.  Fig.  7  gives  the  resulting  stresses 
for  the  first  cycle  and  the  growth  of  the  total  plastic  strain  with  further 
cycling.  During  the  cooling  half  of  the  first  cycle,  tensile  plastic  deforma¬ 
tion  takes  place  at  the  inner  side  of  the  tube  from  x  =  0.7  to  x  =  1,  while 
during  the  heating  half  of  this  cycle  compressive  plastic  deformation  takes 
place  also  at  the  inner  side  from  x  =  0.96  to  x  =  1,  Further  cycling  gives 
exactly  the  same  stress  distributions.  Hence,  the  range  from  x  =  0.96  to 

x  =  1  at  the  inner  side  is  undergoing  cyclic  plasticity  while  no  total  strain 
growth  is  observed.  The  cyclic  steady  state  is  thus  cyclic  plasticity  F^  and 
is  reached  after  the  first  cycle.  Similar  analytical  results  were  obtained 
for  the  same  pressure  and  temperature  difference. 

(e)  For  p  =  900  psi  and  AT  =  50°F.  Fig.  8  gives  the  resulting  stress  dis¬ 
tributions  during  the  first,  tenth  and  twentieth  cycles.  The  total  strain 
growth  is  shown  in  Fig.  12.  It  is  clear,  in  this  case,  that  plastic  deformations 
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occur  on  one  side  only  of  the  tube  during  each  half  cycle,  and  during  the 
subsequent  half  this  side  is  fully  elastic  and  the  other  half  starts  to 
yield.  Such  a  mechanism  is  maintained  during  all  cycles,  thus  causing  a 
growing  total  plastic  strain  over  the  cycles  (ratchetting)  with  a  decreas¬ 
ing  rate  of  growth  until  a  cyclic  steady  state  of  elastic  shakedown  is 
reached.  It  can  be  seen  from  Fig.  8,  that  the  plastic  fronts  y^  and  j2  are 
getting  closer  with  further  cycling  until  they  coincide  at  x  ■  0.56  when 
the  steady  state  is  approached.  This  is  seen  almost  after  18  cycles.  From 
Fig.  12,  it  can  be  seen  that  the  rate  of  growth  of  the  total  strain  is 
large  at  the  beginning  and  is  decreasing  gradually  until  the  steady  state 
is  observed  with  a  value  for  the  total  strain  of  about  0.24Z.  This  is 
considered  the  asymptotic  ratchet  strain  for  the  loads  imposed.  The  observed 
ratchetting  regime  is  R^.  Similar  analytical  results  were  obtained  under  the 
same  conditions. 

(f)  For  p  -  200  psi  and  AT  »  300°F.  The  resulting  stress  distributions 

are  shown  in  Fig.  9  for  the  first  and  almost  all  subsequent  cycles.  The 
total  strain  growth  with  cycling  is  shown  in  Fig.  12.  In  this  case  plastic 
deformations  are  taking  place  on  both  sides  of  the  tube  during  each  he If 
cycle  in  such  a  fashion  as  to  increase  the  total  strain  over  each  half  cycle 
by  a  certain  decreasing  increment  until  a  steady  state  of  reversed  plasticity, 
F^,  is  reached.  For  the  case  considered,  compressive  plastic  deformations 
occur  at  the  outer  side  from  x  *  0  to  x  «  0.38  and  tensile  plastic  deformations 
occur  at  the  inner  tube  side  from  x  *  0.52  to  x  *  1  during  the  cooling  half  of 
any  cycle.  For  the  heating  half  cycle,  the  situation  is  reversed  and  tensile 
plasticity  takes  place  from  x  *  0  to  x  ■  0.48  whilst  compressive  plasticity 
occurs  from  x  *  0  to  x  *  0.66  to  x  *  1.  The  cyclic  steady  state  is  reached 
after  almost  10  cycles  with  an  asymptotic  ratchet  strain  of  -  0.01Z.  The 
observed  ratchetting  regime  for  this  case  is  Similar  analytical  results 

were  found  under  the  same  operating  conditions. 

(g)  For  p  »  900  psi  and  AT  100°F.  The  resulting  stress  distributions  are 
shown  in  Fig.  10  for  the  first,  tenth  and  twentieth  cycles.  The  total  strain 
growth  is  shown  in  Fig.  12.  Here,  a  phase  of  transient  ratchetting  similar  to 
R^  is  seen  up  to  the  seventh  cycle  until  a  sufficient  degree  of  hardening  is 
attained  which  initiated  reversed  plasticity  at  the  inner  tube  side  in  a 


fashion  similar  to  F^.  With  further  cycling  the  plastic  fronts  for  two 
subsequent  half  cycles  are  getting  closer,  and  after  the  heating  half  of 
the  7th  cycle  a  new  compressive  plastic  zone  is  observed  at  the  inner  side 
of  the  tube.  The  size  of  this  zone  is  Increased  with  cycling  till  it 
covers  the  range  from  x  *  0.96  to  x  •  1,  at  the  steady  state.  This  range 
is  therefore  undergoing  cyclic  plasticity  F^.  The  total  accumulated 
plastic  strain  reaches  its  asymptotic  value  of  about  0.65Z  after  25  cycles 
where  the  ratchettlng  plastic  fronts  coincide  at  x  *  0.55.  A  cyclic  steady 
state  of  reversed  plasticity  at  the  inner  tube  side  is  thus  maintained. 

These  results  are  found  to  be  in  a  very  good  agreement  with  analytical 
results  under  similar  conditions. 

(h)  For  p  -  900  psl  and  AT  »  300°F.  The  resulting  stress  distributions  - 
using  the  K.H.  uniaxial  version  -  are  shown  in  Fig.  11  for  the  first  and 
tenth  cycles  (which  was  almost  the  same  for  all  subsequent  cycles).  The 
total  strain  growth  is  shown  in  Fig.  12.  For  this  case,  ratchettlng  is 
observed  at  the  first  few  cycles  in  a  mode  similar  to  until  a  sufficient 
degree  of  hardening  is  attained  to  initiate  reversed  plasticity  at  the  outer 
tube  side  (this  was  observed  after  the  6th  cycle),  thus,  giving  rise  to  a 
ratchettlng  regime  R,  which  will  have  a  steady  state  of  cyclic  plasticity  F£. 
The  cyclic  steady  state  is  reached  almost  after  15  cycles,  with  the  ratchet- 
ting  plastic  fronts  coinciding  at  x  *  0.5  and  reversed  plasticity  covering 
the  ranges  from  x  *  0  to  x  *  0.38  at  the  outer  side  and  from  x  *  0.66  to 
x  *  1  at  the  inner  side.  The  accumulated  total  plastic  strain  at  the  steady 
state  is  about  1.3Z  which  is  to  be  considered  the  asymptotic  value  for  the 
ratchettlng  strain.  The  observed  ratchettlng  regime  is  called  R^.  Analytical 
evaluation  of  stresses,  strains,  ratchettlng  strains  and  asymptotic  strain 
gave  similar  figures  to  those  obtained  numerically. 

4.2.2  Istropic  hardening  results 

(a)  When  the  isotropic  hardening  uniaxial  version  of  the  program  was  used 
for  p  «  100  psl  and  AT  •  50°F,  the  results  were  exactly  equal  to  those 
obtained  from  the  kinematic  hardening  version.  Shakedown  was  observed  after 
the  first  cycle.  The  stress  distributions  and  total  strain  growth  are  there¬ 
fore  as  shown  in  Fig.  4.  The  reason  for  the  similarity  between  kinematic  and 
isotropic  hardening  results  is  due  to  the  fact  that  the  effect  of  the  hardening 
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rule  on  a  virgin  material  can  be  seen  only  when  there  is  more  than  one  half 
cycle  during  which  plasticity  occurs.  Since  in  the  shakedown  regime,  no 
plasticity  takes  place  after  the  first  half  of  the  first  cycle,  it  is  then 
evident  that  the  stresses  and  strains  at  the  end  of  this  half  cycle  are  not 
yet  influenced  by  the  hardening  rule. 

(b)  For  p  ■  100  psi  and  AT  -  100°F.  The  results  obtained  are  also  as 
shown  in *Fig.  5.  The  features  of  the  stress  distributions  over  any  cycle 
and  the  trend  for  cumulative  strain  growth  are  found  to  be  exactly  the  same 
as  those  obtained  from  the  K.H.  version  for  the  same  loads,  though  the  mode 
of  deformation  observed  was  elastic  shakedown  after  the  second  cycle  compared 
with  cyclic  plasticity  for  kinematic  hardening  results.  Ho  plastic  deforma¬ 
tions  were  seen  after  the  second  cycle  since  the  degree  of  hardening  reached 
was  sufficient  to  prevent  the  formation  of  any  new  plastic  strains. 

(c)  For  p  »100  psi  and  AT  «=  200°F.  The  isotropic  hardening  results  are 
coincidentally  similar  to  those  obtained  for  kinematic  hardening  and  a  state 
of  cyclic  plasticity  with  a  slight  decreasing  reversed  plastic  strain  was 
observed  with  further  cycling.  No  total  strain  growth  was  observed,  as  seen 
from  Fig.  6  which  applies  here  also.  Nevertheless  it  is  worth  warning  that 
the  similarity  observed  up  till  this  point  between  the  isotropic  and  kinematic 
hardening  results,  is  not  a  common  rule  for  all  results.  It  is  confined  only 
to  the  stress  distributions  and  trend  of  strain  growth,  but  no  the  regime  of 
deformation,  and  only  for  the  loading  combinations  considered.  For  the  con¬ 
sidered  point  (p  ■  100  psi  ana  AT  ■  200°F),  cyclic  plasticity  at  both  sides 

of  the  tube  was  observed  only  up  to  the  fourth  cycle  after  which  a  state  of 
elastic  shakedown  was  reached. 

From  the  above  three  examples,  it  can  be  generalized  that  for  a  tube 
made  of  an  isotropically  hardening  material  and  for  any  loading  combinations, 
the  cyclic  steady  state  will  always  be  an  elastic  shakedown. 

4.2.3  Results  for  load  conditions  used  in  Corum’s  teat 

For  p  -  700  psi  and  AT  ■  167°F  which  correspond  to  the  conditions  of 
Corum's  test  / 8/,  the  kinematic  hardening  uniaxial  version  of  the  program 
gave  the  stress  distributions  shown  in  Fig.  13  (h  *  0.328)  for  the  first, 
tenth,  twentieth  and  thirtieth  cycles.  The  corresponding  total  strain  growth 
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is  shown  in  Fig.  15.  The  observed  mode  of  deformation  was  ratchetting  E, 
in  which  cyclic  plasticity  was  observed  at  the  inner  tube  side.  It  covers 
the  range  from  x  -  0.9  to  x  -  1  during  the  first  cycle  and  increases 
gradually  till  the  steady  state  is  approached.  The  steady  state  is  almost 
reached  after  25  cycles.  The  ratchetting  plastic  fronts  are  at  x  *  0.48 
and  x  *  0.54  after  30  cycles.  Plasticity  covers  the  range  from  x  »  0.78 
to  x  «  1,  and  the  total  cumulative  strain  is  about  0.73%. 

For  the  same  pressure  (700  psi)  and  temperature  difference  (167°F), 
the  linear  temperature  distribution  (h  *  0)  was  tried  and  the  K.H/uniaxLal 
results  are  as  shown  in  Figs.  14  and  15.  The  mode  of  deformation  observed 
is  also  ratchetting  but  regime  Rj  with  cyclic  plasticity  seen  from  the 
first  cycle  at  both  sides  of  the  tube  covering  the  ranges  from  x  *  0  to 
x  ■  0.02  and  from  x  -  0.98  to  x  *  1.  These  ranges  increase  with  further 
cycling  till  the  steady  state  is  reached  after  almost  20  cycles.  It  covers 
then  the  ranges  from  x  *  0  to  x  *  0.12  at  the  outer  skin  and  from  x  *  0.88 
to  x  *  1  at  the  inner  skin.  The  observed  cyclic  steady  state  is  reversed 
plasticity,  regime  Fj.  The  value  of  the  cumulative  strain  growth  at  the 
end  of  the  30th  cycle  is  0.6%  which  is  lower  than  than  that  for  h  *  0.328. 

When  the  isotropic  hardening  uniaxial  version  was  used  for  Corum's 
test  conditions,  the  stress  distributions  for  the  bilinear  temperature 
distribution  (h  *  0.328)  are  as  shown  in  Fig.  16,  for  the  first,  tenth, 
twentieth  and  thirtieth  cycles.  The  strain  growth  with  cycling  is  shown 
in  Fig.  18.  The  observed  mode  of  deformation  is  ratchetting  with  cyclic 
plasticity  seen  from  the  first  cycle  at  the  inner  tube  side,  covering  the 
range  from  x  *  0.96  to  x  ■  1.  This  range  does  not  very  much  extend  with 
cycling.  After  30  cycles  it  covers  from  x  *  0.94  to  x  *  1.  The  ratchetting 
plastic  fronts  are  clearly  approaching  each  other  with  cycling.  At  the 
first  cycle  they  were  at  x  »  0.46  and  x  *  0.68,  meanwhile  at  the  30th  cycle 
they  were  at  x  *  0.5  and  x  *  0.56.  A  pronouncedly  increasing  ratchetting 
strain  was  still  observed  after  the  30th  cycle  and  no  sign  of  the  steady 
state  was  yet  seen.  The  total  cumulative  plastic  strain  was  1.33%  which  is 
almost  twice  as  much  as  that  obtained  from  kinematic  hardening  results. 

For  the  linear  temperature  distribution  (h  »  0),  the  isotropic  hardening/ 
uniaxial  results  are  as  shown  in  Figs.  17  and  18.  The  deformation  mode  is  pure 
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ratchetting  with  no  reversed  plasticity  zones  at  either  sides  of  the  tube. 

This  ratchetting  node  is  similar  to  The  ratchetting  plastic  fronts 
are  at  x  *  0.36  and  x  «  0.62  at  the  first  cycle  and  generally  approach 
each  other  till  they  reach  x  *  0.5  and  x  *  0.52  after  30  cycles.  The 
steady  state  is  not  clearly  seen  after  30  cycles,  since  the  incremental 
strain  growth  was  still  pronounced.  The  general  trend  of  the  total  strain 
growth  curve.  Fig.  18,  shows  that  the  steady  state  is  expected  after  few 
more  cycles  following  the  30th.  The  total  strain  growth  is  0.8Z,  which 
is  1.25  times  higher  than  the  corresponding  kinematic  hardening  value. 

A  new  set  of  results  for  the  same  test  conditions  is  produced 
using  the  biaxial  tube  program.  This  program  was  tested  for  validity  of 
performance  for  various  loading  points  on  the  Bree  diagram  and  was  found  to 
be  numerically  reliable.  When  the  kinematic  hardening  version  of  this 
program  was  used,  the  results  for  Corum's  test  conditions  were  as  shown  in 
Figs.  19  and  21  for  the  bilinear  temperature  distribution  (h  *  0.328)  and 
in  Figs.  20  and  21  for  the  linear  temperature  distribution.  For  the  first 
case,  the  hoop  stress  (o^)  distributions  during  the  first,  tenth  and 
twentieth  cycles  are  shown  in  Fig.  19,  while  the  growth  of  the  hoop  strain 
at  the  outer  surface  of  the  tube,  e.  (0),  is  as  shown  in  Fig.  21.  Ratchetting 

V 

was  also  observed  here,  with  smaller  ratchetting  strain  increments  and  higher 
rate  of  decrease  of  plastic  strain  increments  with  cycling,  as  compared  to  the 
corresponding  values  (Fig.  15)  obtained  for  the  uniaxial  model  under  similar 
conditions  of  hardening  rule  and  temperature  distribution.  A  cyclic  steady 
state  of  reversed  plasticity  was  observed  after  18  cycles  with  an  asymptotic 
value  of  the  hoop  strain  at  the  outer  tube  surface  of  0.34Z.  This  is  almost 
half  of  that  obtained  from  the  uniaxial  model  results. 

The  linear  distribution  of  temperature  under  the  same  loading  conditions 
(p  *  700  psi  and  AT  *  167°F)  showed  similar  trends.  Ratchetting  was  also 
observed  but  with  lower  values  and  higher  rates  of  the  plastic  strain  increments 
than  those  for  h  *  0.328.  Therefore,  the  steady  state  was  observed  earlier, 
after  15  cycles,  with  an  asymptotic  value  of  the  hoop  strain  at  the 
outer  radius  of  the  tube  of  0.25Z  which  is  almost  3/4  the  corresponding  value 
for  h  ■  0.328.  The  hoop  stress  distributions  during  the  first,  tenth  and 
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twentieth  cycles  are  as  shown  in  Fig.  20,  and  the  hoop  strain  growth  is  as 
shown  in  Fig.  21. 

Similar  to  the  situation  observed  from  the  uniaxial  model  predictions, 
the  isotropic  hardening  predictions,  using  the  biaxial  model,  showed  excessive 
estimation  of  ratchetting  strains,  as  compared  with  kinematic  hardening  predic¬ 
tions  for  both  linear  (h  ■  0)  and  bilinear  (h  *  0.328)  temperature  distributions 
Since  iso'troplc  hardening  assumptions  are  not  in  a  good  correlation  with  actual 
observed  behaviour  of  304SS,  only  the  hoop  strains  at  the  outer  tube  surface  - 
being  of  prime  importance  -  are  shown  in  Fig.  21.  The  distribution  of  hoop 
stresses  during  certain  cycles  is  similar  in  form,  but  different  in  magnitude, 
to  those  shown  in  Figs.  19  and  20. 

As  could  be  seen  from  Fig.  21,  the  total  cumulative  hoop  strain  at  the 
outer  tube  surface  after  30  cycles,  is  0.62  for  h  ■  0.328,  and  0.322  for  h  *  0. 
It  is  seen  also  tbat  the  steady  state  was  not  reached  for  h  *  0.328  and  was 
reached  after  almost  25  cycles  for  h  »  0.  A  summary  of  all  the  test  point 
results  are  given  in  Table  III. 

For  comparison  with  experimental  results,  the  most  suitable  sets  of 
numerical  results  would  be  those  obtained  when  kinematic  hardening  assumptions 
together  with  the  bilinear  temperature  distribution  were  adopted.  This  is 
simply  because  the  hardening  characteristics  of  type  304  stainless  steel  are 
better  approximated  by  kinematic  hardening  rule  than  by  isotropic  hardening, 
and  the  actual  thermal  gradient  is  better  approximated  by  a  bilinear  distribu¬ 
tion  (with  h  «  0.328)  rather  than  a  linear  one.  Fig.  22  shows  the  experimental 
and  numerical  results  for  13  complete  cycles.  The  numerical  predictions  are 
those  obtained  from  both  the  uniaxial  and  biaxial  models  of  the  tube,  under 
kinematic  hardening  assumptions  and  a  bilinear  temperature  distribution.  It 
is  clear  that  the  uniaxial  model  predictions  excessively  overestimate  the 
ratchet  strains,  whilst  the  biaxial  model  predictions  slightly  underestimate 
then.  This  argument  applies  also  to  the  plastic  strain  increment  over  each 
cycle.  After  13  cycles,  the  uniaxial  model  gives  a  cumulative  strain  of  1.95 
times  the  actual  observed  strain  and  a  plastic  strain  Increment  of  3  times  the 
experimental  value,  meanwhile  the  biaxial  model  gives  a  cumulative  strain  of 
only  0.89  times  the  experimental  and  a  plastic  strain  increment  of  only  1.22 


times  the  experimentally  observed  plastic  strain  increment.  It  is  therefore 
evident  that  the  biaxial  model  gives  much  better  predictions  for  the  tube 
problem  than  the  uniaxial.  The  small  differences  between  the  experimental 
results  and  the  available  biaxial  model  predictions  may  be  due  to  the  combined 
effect  of  the  bilinear  approximation  of  the  actual  thermal  gradient,  the 
elastic-linearly  plastic  assumption  as  the  flow  rule  for  the  tube  material, 
the  kinematic  hardening  assumption  as  the  hardening  rule  for  the  tube  material, 
and  the  neglection  of  changes  in  material  properties  when  its  temperature  is 
changed.  It  is  expected  that  better  predictions  will  be  obtained,  if  the  above 
considerations  are  more  accurately  Implemented  in  our  numerical  procedure. 

SECTION  5 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  problem  of  a  tube  subjected  to  steady  internal  pressure  and  cyclic 
radial  thermal  gradient  has  been  investigated.  A  new  numerical  procedure  was 
constructed  to  solve  the  same  problem.  The  validity  of  this  numerical  proced¬ 
ure  was  extensively  verified  through  comparisons  with  available  analytical 
solutions,  namely,  that  based  on  the  kinematic  hardening  assumption  for  an 
elastic-linearily  plastic  tube  material  and  the  uniaxial  model  of  the  tube. 

The  numerical  technique  has  been  further  elaborated  to  accommodate  the 
isotropic  hardening  rule  Instead  of  the  kinematic  hardening  and  to  utilize  a 
more  generalized  biaxial  tube  model  Instead  of  the  uniaxial  one.  The  main 
results  of  this  work  are  summarized  again  hereafter  as  main  conclusions. 


(a)  The  adopted  numerical  procedure  is  a  successful  one.  Its  predictions, 
in  the  case  of  uniaxial  model  with  kinematic  hardening  assumptions, 
were  in  a  very  good  agreement  with  analytical  results.  Expected 
regimes  of  deformation  were  exactly  as  predicted  from  the  full  solution 
map  of  Fig.  3.  Asymptotic  ratchet  strains  were  also  the  same  as  those 
obtained  from  the  same  figure.  The  stresses  at  any  point  of  the  tube 
thickness  during  any  cycle  and  the  values  of  the  plastic  fronts  were 
almost  the  same  as  calculated  analytically. 


(b)  The  numerical  solution  Is  not  confined  to  the  kinematic  hardening 
assumption  only,  hut  also  isotropic  hardening  was  made  to  replace 
the  kinematic  one,  for  which  no  analytical  solutions  are  available, 
and  it  was  possible,  therefore,  to  reach  the  thermo-elasto-plastic 
solution  for  the  tube  problem  under  any  loading  conditions.  Shake- 
down,  reversed  plasticity  and  ratchetting  were  also  observed  here 
with  a  common  cyclic  steady  state  of  shakedown  for  all  regimes. 

This  is  due  to  the  fact  that  under  isotropic  hardening  assumption 
the  size  of  the  yield  surface  is  increased  with  the  increasing 
number  of  load  cycles  until  a  sufficient  degree  of  hardening  is 
attained,  such  that  the  loads  are  no  longer  able  to  cause  plasticity 
and  fully  elastic  behaviour  prevails  then  (shakedown) . 

(c)  A  new  model  has  been  developed,  for  which  the  actual  stress  and 

strain  conditions  were  maintained.  This  model  has  been  found  more 

accurate  and  superior  to  the  uniaxial  model  initially  proposed  and 

that  has  been  widely  adopted  in  most  of  the  trials  to  solve  the  tube 

problem.  The  biaxial  tube  model  is  very  close  to  actual  conditions 

where  all  stresses  and  strains  are  present  in  the  three  dimensions, 

i.e.  a  ,  a„t  a  ,  a  ,  en  and  e  .  Since  the  tube  problem  is  one 
r  o  z  r  0  z 

in  which  axial  symmetry  exists,  all  stresses  and  strains  -  except  e  ~ 

z 

were  dependent  on  the  radius  only.  The  axial  strain  is  constant, 
since  the  tube  was  assumed  to  be  in  a  state  of  generalized  plane 
strain.  Due  to  its  faithful  representation  of  the  tube,  the  results 
obtained  using  the  biaxial  model  has  been  found,  as  expected,  much 
closer  to  actual  experimental  values.  Also,  the  biaxial  model  predic¬ 
tions  of  ratchet  strains,  when  isotropic  hardening  was  assumed,  were 
well  above  those  obtained  from  the  same  model  but  with  kinematic  hard¬ 
ening  assumption.  Fig.  21;  a  situation  similar  to  that  observed  for 
the  uniaxial  model. 

(d)  The  effect  of  the  hardening  rule  on  the  ratchet  strain  and  number  of 
cycles  needed  to  reach  the  steady  state  has  been  found  so  pronounced. 
Isotropic  hardening  assumption  has  resulted  in  greater  ratchet  strains 
than  those  produced  under  the  kinematic  hardening  assumption  for  both 
the  uniaxial  and  biaxial  tube  models.  Also  the  tendency  to  reach  the 
cyclic  state  has  been  found  to  be  at  a  relatively  smaller  number  of 


thermal  cycles  under  kinematic  hardening  rule  than  those  for  the 
isotropic  one.  This  vas  true  for  the  uniaxial  as  well  as  the 
biaxial  tube  model. 

(e)  The  effect  of  'h',  the  bilinear  intersection,  on  the  ratchet  strains 
and  the  arrival  at  the  steady  state  was  also  a  very  clear  one.  For 
both  isotropic  and  kinematic  hardening  and  for  both  the  uni.axial  and 
biaxial  tube  models  the  ratchet  strains  were  relatively  higher  in 
the  case  of  a  bilinear  temperature  distribution  (h“0.328)  than  that 
for  the  linear  distribution  (h  *  0).  The  cyclic  steady  state  was 
arrived  at,  considerably  earlier  in  the  case  of  a  linear  distribution 
than  a  bilinear  one.  Table  III  shows  the  above  effect  quantitatively. 

(f)  The  difference  found  between  the  experimental  and  biaxial  model 
results  (the  latter  being  about  10Z  lower  than  the  former),  is  not 

as  dreadful  as  those  observed  when  Bree’s  uniaxial  tube  model  was  used 
(about  95Z  above  the  actual  results).  Possible  sources  of  error  in 
the  load  and  material  modeling  are  : 

i.  the  elastic-linearly  plastic  assumption  with  6  =  0.0253  for  the 
tube  material, 

ii.  the  kinematic  hardening  assumption  for  the  304SS, 

iii.  the  bilinear  temperature  distribution  with  h  *  0.328,  which  is 
only  an  approximation, 

iv.  the  constancy  of  material  properties  assumption,  which  neglects 
the  effect  of  temperature  on  properties. 

Nevertheless  it  is  worth  mentioning  that  the  adopted  numerical  technique 
can  handle  ary  type  of  flow  rule,  hardening  rule,  temperature  distribution  and 
any  set  of  data  describing  the  change  of  material  properties  with  temperature, 
once  these  are  available.  Even  the  control  history  can  be  changed  from  cycle 
to  cycle.  Also,  any  tube  dimensions  can  be  used,  including  thick-walled  tubes, 
since  the  biaxial  model  simulates  basically  thick-walled  tubes.  The  full 
solution  map  for  a  certain  tube  (Bree  diagram)  can  then  be  numerically  construe 
ted  for  any  set  of  material  data. 


From  the  foregoing,  it  Is  clear  that  the  flow  rule  and  the  hardening 
rule  play  a  very  important  role  in  the  deformatlonal  pattern  of  the  struc¬ 
ture.  Although  the  kinematic  hardening  rule  reproduces  the  observed 
Bauschinger  effect  upon  stress  reversals,  the  cyclic  steady  state  under 
conditions  of  fully  reversed  strains,  is  reached  during  the  first  cycle. 

On  the  other  hand,  isotropic  hardening  theory  reaches  a  purely  elastic 
cyclic  state  after  a  finite  number  of  cycles  and  does  not  show  a  Bauschinger 
effect.  The  two  theories  are  seen  to  provide  two  extreme  limits  of  behaviour, 
between  which  actual  material  response  as  well  as  all  remaining  hardening 
theories  lie.  Fig.  23.  In  this  respect,  various  efforts  have  been  done 
recently  towards  a  more  accurate  description  of  material  behaviour.  Among 
the  most  important  models  presented  in  the  literature  are  those  of  Mroz  / 10/ , 
Eisenberg  /ll/  and  Daf alias  and  Popov  / 12/ • 

Mroz  / 10/  proposed  a  material  model  which  utilizes  a  field  of  work¬ 
hardening  moduli  as  shown  in  Fig.  24.  The  monotonic  stress-strain  curve  is 
approximated  by  a  number  of  linear  segments.  During  monotonic  loading,  the 
stress  increases  in  steps  and  plastic  strain  is  defined  by  the  appropriate 
work-hardening  modulus. 

Eisenberg  / 11/  assumed  the  universal  applicability  of  Ramberg-Osgood 
relations,  e  =  a/E  +  (a/k)n  for  the  monotonic  and  subsequent  branches  of  the 
hysteresis  loops  and  developed  growth  laws  for  the  material  parameters  k  and 
n  on  the  basis  of  observed  hardening  of  304SS. 

Daf alias  and  Popov  / 12/  introduced  the  concept  of  a  limit  or  bounding 
surface  which  encloses  the  yield  surface  in  stress  space,  Fig.  25.  The  limit 
surface  can  translate,  shrink  or  expand  but  never  intersect  with  the  yield 
surface.  The  plastic  tangent  depends  on  the  proximity  between  the  two  surfaces 
at  the  point  of  loading.  The  material  exhibits  a  memory  of  past  loading  history 
by  means  of  plastic  internal  variable  and  a  projected  forsight  by  means  of  the 
proximity  between  actual  and  image  stress  points.' 

It  is  recommended  that  these  and  other  models  should  be  used  in  further 
work  related  to  the  tube  problem. 
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TABL^  I  -  Internal  pressure  -  temperature  difference  combinations  and 

corresponding  Increment  of  t  emperature  difference  and  expected 
mode  of  behavior  used  for  in  specting  the  uniaxial  program  with 
kinematic  hardening  assumption 


P  (psi) 

100 

100 

100 

— 

100 

900 

200 

900 

900 

AT  (°F) 

50 

200 

300 

ICO 

50 

300 

100 

300 

AT'(°?) 

20 

20 

30 

5 

20 

30 

20 

30 

moce  * 

F1 

F2 

F:i 

R2 

R3 

R< 

TABLE  II  -  Internal  pressure  -  tempera  .ure  difference  combinations  and 
corresponding  increments  of  temperature  difference  used  in 
the  uniaxial  program  with  inotropic  hardening  assumption 


P 

100 

100 

100 

7C0 

AT 

50 

100 

200 

167 

AT' 

o 

CM 

20 

20 

10 

TABLE  Ill-Summary  of  numerical  results  for  the  load  conditions 
(p  *  700  psi  and  AT  ■  167CF)  used  by  Corum  in  his 
Cube  ratchetting  tests 


Model  Uniaxial  Model  Biaxial  Model 


Hardening  rule  K.H.  I,H.  I  K.H.  I.H. 


Fig.  ?  Simplified  flow  chart  for  the  biaxial  program 


Fig.  5  Cyclic  stress  distributions  and  strain  growth  for  P  =  100  psi 
and  At  =  200°F  (Regime  F^) 


Fig.  6  Cyclic  stress  distributions  and  strain  growth  ;or  p  =  100  psi 
and  AT  =  300°F  (Regime  F2) 


Cycles 


Fig.  7  Cyclic  stress  distribution  and  strain  growth  for  P  =  100  psi 
and  T  =  100°F  (Regime  Fj 


n»l 


Fig.  9  Cyclic  stress  distributions 
for  p  =  200  psi  and  AT  =  30 
(Regime  R  ) 
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Fig.  8  Cyclic  stress  distributions  for 
p  =  900  psi  and  At  =  50°F 
(Regime  R.) 
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Fig.  10  Cyclic  stress  distributions  for 
p  *  900  psi  and  AT  *  100  F 
(Regime  R3) 


Fig.  11  Cyclic  stress  distributions 

for  p  -  900  psi  and  AT  -  300°* 
(Regime  R^) 
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Fig.  13  Cyclic  stress  distributions  for  p  »  700  psi 
and  AT  »  167°F  (K.H.  -  uniaxial  version  with 
h  «  0.328) 


V-*. 


Fig.  20  Cyclic  hoop  stress 

distribution  for  700  psi/ 
167°F  (K.H.  -  biaxial 
version  with  h  *  0) 


Rotchet  strain 


Fig.  21  Ratchet  hoop  strain  growth  for  700  psl/167°F  for  both  K.H 
and  I.H.  biaxial  programs  with  h  *  0  and  0.328 


Fig.  22  Uniaxial,  biaxial  and  experimental  results  for  p  *  700  psi 
and  AT  «  167°F 


